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FOREWORD 

This  research  deals  with  the  mechanism,  properties  and 
characteristics  of  vortex  shedding  vibrations  of  bridge  members 
for  the  purpose  of  disseminating  the  researched  information  into 
a  design  aid  in  order  to  mitigate  the  damages  done  to  long,  slender 
members  through  fatigue  caused  by  wind-induced  vibrations. 

The  main  results  of  this  investigation  are  reduced  into 
twelve  sections,  which  are  listed  as  the  texts  of  this  report; 
complicated  formulas  that  require  more  than  ordinary  calculation 
efforts  are  presented  in  graph  form  in  Appendix  A.   The  text 
and  Appendix  are  integral  parts  of  a  proposed  design  guide.   To- 
gether ,  they  can  be  used ,  on  our  opinion ,  by  structural  designers 
in  a  routine  manner. 

To  provide  adequate  background  of  this  design  aid  four 
supplements  are  appended  with  the  text.   These  supplements 
should  be  consulted  by  those  who  wish  to  understand  the  basis 
upon  which  the  design  aid  has  been  developed. 

Supplement  Part  A  is  a  digest  of  literature  which  was  re- 
viewed during  the  course  of  this  investigation.   It  provides 
comprehensive  background  material  on  the  historical  development, 
origin  and  behavior  of  vortex-shedding  mechanism.   These  form  a 
basis  for  establishing  the  methods  and  devices  to  suppress  the 
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formation  of  the  vortex  street  and  thus  reduced  the  severity 
of  vibration  of  the  structural  members.   The  available  Strouhal 
numbers  of  the  most  commonly  known  cross-shapes  of  structural 
applicability  are  listed  as  are  other  useful  engineering  data 
including  damping,  lift  and  end-fixity  coefficients. 

Supplement  Parts  B  and  C  consist  of  all  mathematical 
theories  and  implementation  of  these  theories.   For  the  convenience 
of  the  users ,  we  have  segregated  the  mathematical  and  computa- 
tional aspects  from  each  other. 

The  mathematical  theory  of  free  vibration  including  the 
frequency  equations  and  mode  shapes  are  given  in  Chapter  I , 
Part  B.   The  flexural  vibration  theory  is  applicable  to  a  beam 
with  elastic  constrained  ends  and  with  axial  tension.   This  is 
believed  to  be  the  most  general  formulation  in  the  framework 
of  Euler-Bernoulli  theory.   The  torsional  vibration  theory  is 
applicable  to  thin  closed  and  open  sections  with  consideration 
to  warping  and  axial  tension.   Neither  theory  has  appeared  in 
the  literature  in  this  generality  and  this  is  believed  to  have 
advanced  the  state  of  the  art.  The  detailed  evaluation  of 
natural  frequencies  of  various  types  of  vibration  is  given  in 
Part  B,  Chapter  I.   All  results  are  given  in  graph  form,  in 
Appendix  A  of  the  main  text  and  also  in  Part  E. 
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Numerical  examples  which  demonstrate  how  the  theory  and 
natural  frequency  graphs  can  be  used  to  determine  the  corres- 
ponding wind  velocities  which  cause  resonance  in  the  bridge  mem- 
bers are  given  in  Part  C,  Chapter  I.   The  mode  shapes  of  a  vibrating 
beam  are  used  for  random  vibration  as  described  later  in  the  text. 

The  mathematical  theory  which  determines  the  amplitude  of 
vibration  of  a  rigid  beam  in  single  frequency  resonance  (sub- 
critical  flow) is  given  in  Part  B,  Chapter  II.   The  forcing  func- 
tion was  assumed  to  be  a  van  der  Pol  oscillator,  used  by  Hartlen 
and  Currie  for  a  circular  cylindrical  bar.   In  the  present  work, 
we  adapted  it  to  an  H- section  and  calibrated  the  formulation  to 
the  experimental  results  by  Narita.   The  implementation  of  this 
theory  to  an  elastic  beam  is  given  in  Part  C,  Chapter  II, 
accompanied  by  a  numerical  example. 

The  response  of  a  beam  to  random  vibrations  is  derived  in 
Part  B,  Chapter  III,  in  a  self-content  manner.   The  theory  is 
implemented  to  a  beam  with  constrained  ends.   The  implementa- 
tion of  this  theory  is  given  in  Part  C,  Chapter  III.   The  mean  de- 
formation and  the  mean  squared  deformation  are  obtained  from 
the  formulations  given  in  Part  B,  Chapter  I.   The  recommended 
lift  coefficients  are  suitable  for  H-beams  and  box  beams.   It 
was  shown  that  in  most  cases  the  one  term  solution  will  give 
better  than  95%  accuracy.   Since  the  first  mode  of  vibration  is 


-4- 


consistent  in  sign,  the  rms  deflection  is  the  same  as  the  mean 
deflection.   For  this  reason,  the  flexural  stress  as  shown  in 
Part  C,  Chapter  II,  can  be  determined  by  a  simpler  procedure 
than  the  conventional  double  differentiation  method  which  is 
quite  complicated  in  computation. 

Supplement  Part  D  consists  of  a  comprehensive  bibliography 
of  nearly  900  entries ,  accompanied  by  a  summary  of  search  proce- 
dure and  leading  articles.   Supplement  Part  E  consists  of  a 
complete  set  of  graphs  for  natural  frequency  determination 
especially  for  unsymmetrical  cases,  i.e.,  beams  having  unequal 
end- fixities. 
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1.   Notations 

A  =  Cross-sectional  area 

B  =  Amplitude  of  nth  cycle  of  vibration 

CL  =  Lift  Coefficient 

E  =  Young's  modulus  of  elasticity  of  steel 

=  29  x  106  psi 

I  =  Area  moment  of  inertia  of  member  under  consideration 

1^  =  Area  moment  of  inertia  of  lateral  members  in  same  joint 

L  =  Length  of  member  under  consideration 

Lj  =  Length  of  member  in  same  joint 

Re  =  Reynolds  number 

S_  =  Power  spectrum  function 

St  =  Strouhal  number 

V  =  Wind  velocity 

a  =  Moment  coefficient  =  M/qX 

b  =  Deflection  coefficient  =  w   EI/qL^ 

max 

c  =  Extreme  fiber  distance 

d  =  Depth  of  member  in  direction  perpendicular  to  wind 
direction 

e  =  Naperian  base  of  logarithm 

f j  =  ith  natural  frequency 

f  =  Frequency  in  Hertz 

f^  =  Flexural  (bending)  stress 

h  =  width  of  member 

Pa  =  mass  density  of  air 

f)^  =  mass  density  of  member 
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q  =  Uniformly  distributed  transverse  load 

w  =  Deflection 

Wmax  Maxumum  deflection  (amplitude) 

w  =  Mean  deflection  assumed  to  be  h   w 

max 

Y     =  Structural  damping  coefficient 

=  <S/2tt 

6  =  Logarithmic  decrement 

a     =   End- fixity  condition 

w  =  Circular  natural  frequency 
2.   Scope 

This  chapter  covers  the  design  of  slender  bridge  members 
subject  to  wind-induced  vibrations  characterized  by  excursions 
in  the  direction  perpendicular  to  the  wind  direction.   The 
vibrations  considered  herein  are  caused  by  vortex-shedding  in 
both  periodical  (sub-critical)  and  random  (super-critical)  ranges 
of  Reynolds  numbers.   Vibrations  due  to  flutter,  buffeting,  and 
galloping  are  beyond  the  scope  of  this  chapter.   The  structural 
members  are  restricted  to  those  of  which  the  cross  section 
has  at  least  two  axes  of  symmetry,  e.g.,  circular,  rectangular, 
box,  H-sections,  double -channel  sections,  etc.   These  type  of 
sections  are  commonly  used  as  tension  members  in  the  modern 
bridges.   Compression  members  employing  built-up  sections  usually 
have  only  one  axis  of  symmetry  and  are  capable  of  coupled 
torsional  and  flexural  vibrational  nodes.   A  comprehensive  analysis 
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must  be  conducted  specifically  for  such  cases. 

3.  Assumptions 

(a)  The  assumptions  made  in  a  general  "strength  of  material" 
analysis  of  a  member  subjected  bending  and  torsional 
conditions  shall  apply  to  the  member  being  analyzed. 

(b)  Wind  forces  are  assumed  to  be  acting  in  the  lateral 
direction  of  the  member  and  perpendicular  to  the 
principal  axes  of  the  member. 

(c)  Members  are  assumed  to  be  spaced  such  that  there  is 
sufficient  room  for  the  wake  to  fully  develop  without 
interference  from  any  adjacent  members,  behind  the 
member  being  analyzed. 

(d)  It  is  assumed  that:   (1)   The  natural  frequencies  of  the 
member  being  analyzed  are  widely  separated  from  one 
another,   (2)   The  damping  is  small,  and  (3)  The  lift 
forces  are  spatially  independent. 

4.  Strouhal  Numbers 

For  design  purpose,  the  Strouhal  numbers  for  common  cross- 
sections  shall  be: 


Circle 

Rectangle 
h 

H-Section 

0.20 

0.13(_<1.0) 
d 

0.12 
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5.   Damping 

Damping  coefficient  shall  be  expressed  in  terms  of  the 
logarithmic  decrement,  which  is  defined  to  be  the  natural 
logarithm  of  the  ratio  of  two  successive  peak  amplitudes  in 
a  free,  decreasing  oscillation,  i.e.,  if  Bn  and  Bn+^  are  the 
nth  and  (n+l)th  amplitudes,   respectively,  the 

Bn 
6  =  Log 


e  „ 
Bn+1 

For  design  purposes,  6  may  be  taken  as  0.02  for  members  capable 
of  resisting  flexural  moments, or  0.08  for  flexible  members 
including  steel  wire  ropes.   Structural  damping  may  be  taken  to 

be  y  =  6/27T 

6.   Lift  Coefficient 

The  fluctuating  or  time -dependent  lift  coefficient  is 
dependent  on  the  Reynolds  number  and  surface  characteristics  of 
the  member.   For  accurate  analysis  the  lift  coefficient  should 
be  determined  by  wind  tunnel  tests.   For  design  purposes, 
0.6  and  1.2  may  be  used  as  the  r.m.s.  and  peak  values  of  lift 
coefficients,  respectively. 

7.   End-Fixity  Condition  of  a  Member 

The  end  fixity  condition  of  a  long  slender  member  in  a  truss  or 
a  frame  is  given  by 

l    I   I./L. 

H    ill 

o  =   -      


n     I/L 
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Where  I   Ii/Li  =  Sum  of  stiffness  of  all  the  lateral  members 
i 

in  a  joint  other  than  member  under  consideration 

I/L  =  Stiffness  of  the  member  under  consideration 

n  =  A  constant  indicating  the  continuity  condition  at 

the  joint,  i.e.,  n  =  1  for  members  which  terminate  themselves  at  the 

joint,  such  as  a  hanger;  n  =  2  for  members  which  are  continuous  through 

the  joint  such  as  the  lower  chord  of  a  truss. 

8.  Natural  Frequency 

The  natural  frequencies  of  a  bridge  member  is  affected  by 
the  section  properties,  axial  tensile  force* and  the  end-fixity 
of  the  member.   The  end-fixity  of  a  member  shall  be  determined  by 
the  method  as  described  in  Section  1_     of  these  recommendations.   The 
recommended  prevailing  tensile  force  in  the  member  shall  be 
that  caused  by  combination  of  dead  load  and  one  fourth  of  the 
design  live  load.   The  natural  frequency  coefficients  for  the 
first  three  f lexural  and  torsional  modes  may  be  obtained  from 
Figs.  1  through  10,  for  all  cases  except  for  thin  closed  sections. 
In  the  latter  case,  natural  frequency  coefficients  shall  be 
obtained  from  Table  1. 

9.  Nominal  Wind  Velocity  at  Resonance 

The  nominal  wind  velocity  which  causes  resonance  in  a 

bridge  member  shall  be  determined  by 

*The  figures  should  be  used  for  tensile  members  only.   The  general 
method  is  applicable  to  compression  members  for  which  the  natural 
frequencies  are  lower. 
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fd 
v  = 


St 

This  velocity  shall  be  used  to  compute  the  Reynolds  numbers 

associated  with  the  flow.   Reynolds  number  may  be  computed  by 

R  =  780  Vd 
e 

if  V  is  in  MPH  and 
d  in  inches. 

The  Reynolds  number  will  indicate  the  flow  regime  for  circular  members, 
If  the  Reynolds  number  is  below  1.0  x  10  ,  the  flow  is  in  the  Strouhal 
range.   If  the  Reynolds  number  is  within  1.0  x  10J  to  4.0  x  10   the 
flow  is  in  the  transition  range  depending  on  upstream  turbulence.   If 
the  Reynolds  number  is  above  4.0  x  10  ,  the  flow  is  in  the  super- 
critical range  and  random  shedding  takes  place.   This  information 
shall  be  used  to  determine  which  of  the  two  approaches  described  in 
the  next  two  sections  shall  be  used  to  calculate  the  deflection  in 
the  member  For  structural  shapes,  subcritical  flow  prevails  for  all 
Reynolds  numbers. 

10.   Deflection  in  Subcritical  Flow 

The  deflection  of  the  bridge  member  in  subcritical  flow 

shall  be  analyzed  by  Fig.  11  which  shows  a  relation  between 

the  non-dimensional  amplitude  w/d  and  a  non-dimension  damping 
YAP* 


Thus,  the  mean  deflection  w  can  be  determined  from  known 
Pad2 

area,  depth  and  structural  damping  coefficient.  The  corresponding 
flexural  stress  can  be  computed  by  a  method  suggested  in  section  12. 
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11.   Deflection  of  Member  in  Supercritical  Flow  (Random  Shedding) 

Condition 

The  mean  deflection  w  of  a  bridge  member  undergoing  a 
random  vibration  shall  be  determined  by 


w  _6      V      /  Sf  (Stj)  Si^TT 

-=  3.1  xio"   CL 

d  A(f,)2/       6 


Where  S   (St) is  the  normalized  power  spectrum.   Power  spectra  for  cylindrical, 
square  and  H-sections  are  given  in  figs. 12,  13A  and  13B,  respectively. 
The  formula  is  valid  provided  that  fi/f^  is  less  than  0.2.   The 
corresponding  flexural  stress  can  be  determined  by  a  method 
suggested  in  the  next  section. 

12.   Bending  Stress  in  Member  Undergoing  a  Known  Deflection 
The  flexural  stress  caused  by  wind  induced  vibrational 

deformation  shall  be  calculated  by: 

f  _  a   w  d  2 
fb  ~  5-  E  a"  <l  } 

where  a/b  is  given  by  Fig.  14.  corresponding  to  a  known  end 
fixity  condition  c,  as  determined  in  Section  7. 

This  bending  stress  shall  be  added  to  the  stresses  imposed 
on  the  member  by  the  normal  dead  and  25  percent  of  live  load 
combination.   If  the  total  computed  stress  exceeds  the  endurance 
limit  of  the  structural  material,  the  designer  should  employ  a 
suitable  measure  to  obviate  this  condition. 
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APPENDIX  A 
The  following  notations  are  given  for  Figs  1-10. 

For  bending  vibrations  (open  and  closed  sections) 

2 
TL 

p  =  j — 


EI 
where    T  =  Axial  tensile  force  in  the  member 
L  =  Span  length 

E  =  Young's  modulus  of  elasticity 
I  =  Moment  of  inertia  (transverse  to  wind  direction) 

pAco  L4 

Q  =  m-^ 

EI 

where    p  =  Mass  density  of  member 
m 

A  =  Cross-sectional  area  of  member 
o)  =  nL11  circular  natural  frequency 
For  torsional  vibration  (thin  open  sections) : 


T      2 
(GC  +  -   I  )L 

P  =  — - 


A   _-P_ 


EC 
w 


2   4 
P  I  to   L 


E   C 
w 

where     G  =  Modulus  of  elasticity  for  shear 
C  =  Torsional  constant 


I  =  Polar  moment  of  inertia 
P 

C  =  Warping  constant 
For  torsion  of  thin  closed  sections,  see  Table  1. 
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Fig,  1  Natural  frequency  at  Different  Tension 
Symmetrical  Beam  with  Hingbd  Ends  (<j=0) 
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Fig,  3     Natural  Frequency  at  Different  Tension 
Symmetrical  Beam  with  <?=10 
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Fig.  4  Natural  Frequency  at  Different  Tension 
Symmetrical  Beam  with  ct=5G 
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Fig.  5     Natural  Frequency  at  Different  Tension 
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Fig.  6  Natural  Frequency  at  Different  Tension 
Symmetrical  Beam  with  a =500 
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Fig.  7  Natural  Frequency  at  Different  Tension 
Symmetrical  Beam  with  <j=] 
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Fig.  8  Natural  Frequency  at  Different  Tension 
Symmetrical  Beam  with  a  =10,000 
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Symmetrical  Beams  with  Fixed  Ends  (<*-*») 


-24- 


Table  1.  Natural  frequency  coefficients  for  thin  closed  sections 


1 

Mode  No. 

Natural  Frequency  Coefficients ,  ^/Q 
For  Any  P 

1 
2 
3 

1.77 
2.51 
3.07 
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PART  A 
GENERAL  DISCUSSIONS 

I .   INTRODUCTION 


In  the  current  design  practice ,  a  dynamic  wind  load  is 
accounted  for  by  converting  it  to  its  static  equivalent.   The 
effect  of  gusts  or  even  hurricane  forces  can  be  safely  taken  care 
of  in  this  manner.   The  wind-induced  vibrations  of  the  bridge 
decks,  especially  the  stiffening  trussess  in  suspension  bridges, 
have  also  attracted  considerable  attention  and  are  carefully 
suppressed  by  stiffening  the  structure,  and  by  modification  of 
cross-sections  with  streamline  fairings  and  other  devices.   One 
problem  which  has  so  far  been  neglected  is  the  wind-induced 
vibration  of  individual  members.   The  individual  members  can  be 
set  in  vibration  by  wind  usually  through  the  vortex  shedding 
phenomenon.   If  the  exciting  frequencies  happen  to  coincide  with* 
natural  frequencies  of  the  structural  members,  rather  severe 
vibrations  can  take  place,  even  in  moderate  wind  velocities  and 
intensities.   Ironically,  it  is  often  the  moderate  wind  that  can 
cause  more  damages  because  of  its  more  commonplace  occurrence, 
rather  than  the  severe  ones,  which  are  rare  events.   For  instance, 
it  is  conceivable  that  a  mild  wind  of  20  knots  can  occur  perhaps 
100  days  per  year  ..   A  typical    bridge  member  in  question 
would  undergo  stress  reversals  of  nearly  2  x  108  cycles  per  year! 
Fatigue  cracks  can  easily  develop  due  to  this  within  a  short 
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period  of  time,  especially  near  the  holes,  re-entrant  corners  and 
other  areas  of  stress  concentrations.   Violent  wind  induced 
vibrations  have  been  actually  observed  in  many  bridge  members  and 
occurrences  of  these  cracks  of  this  type  have  been  discovered 
in  a  number  of  them.   This  report  attempts  to  analyze  the  causes 
and  summarize  possible  remedies  to  this  problem. 
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II.   HISTORICAL  BACKGROUND 

The  fact  that  a  taut  wire  can  be  induced  into  vibration 
by  a  wind  stream  was  known  since  antiquity,  as  evidenced  by  the  use 
of  the  aeolian  harp  as  a  music  instrument  by  the  Greeks  as  early 
as  300  B.C.   The  general  phenomenon  was  first  investigated  by 
Strouhal  in  1878.   In  his  experiments  a  taut  wire  was  attached 
to  a  revolving  frame  about  a  parallel  axis.   It  was  found  that 
the  pitch  of  the  aeolian  tone  was  independent  of  the  length  and 
of  the  tension  of  the  wire,  but  to  vary  with  its  diameter  and 
the  speed  of  the  motion  relative  to  the  air.   He  also  observed 
that  the  intensity  of  sound  was  greatly  enhanced  when  the  speed 
was  such  that  the  aeolian  tone  coincided  with  one  of  the  "Proper 
tone"  of  the  wire  and  that,  other  conditions  being  the  same,  the 
rise  of  temperature  in  the  air  caused  a  slight  but  appreciable 
fall  in  pitch. 

It  had  been  generally  assumed  that  the  excursion  of  the 
wire  would  be  in  the  direction  of  the  wind  stream.   It  was  Lord 
Rayleigh  who  placed  a  taut  wire  near  the  exit  of  chimney  draught 
and  observed  that  the  vibrations  were  effected  in  a  plane  per- 
pendicular to  the  direction  of  the  wind,   The  existence  of  a  vor- 
tex sheet  was  known  to  Leonardo  da  Vinci  in  th.e  fifteenth  century 
who  sketched  a  row  of  vortices  in  the  wake  of  a  bluff  body  in  a 
stream.   Lord  Rayleigh  associated  the  production  of  the  tone  with 
the  vortex  sheets.   He  also  explained  that  the  effect  of  tempera- 
ture on  the  aeolian  tone  was  due  to  the  change  of  viscosity  (which 
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implied  that  the  Strouhal  frequency  must  be  a  function  of  the 
Reynolds  number;  of  course,  the  latter  term  was  unknown  to  him). 

The  wake  behind  a  bluff  body  was  analyzed  by  Helmholtz 
and  Kirchhoff,  but  the  results  were  not  in  agreement  with  the  ex- 
perimental findings  by  Prandtl.   An  independent  approach  was 
used  by  Karman  and  Foppl,  who  left  the  origin  of  the  vortex 
sheet  unanswered  but  addressed  to  the  question:  what  are  the 
stable  arrangements  of  vortices  behind  a  bluff  body?  Karman 
found  that,  for  infinitely  long  straight  vortex  streets  in  two 
rows ,  the  only  stable  arrangement  is  when  the  ratio  between  the 
row  spacing  and  the  spacing  between  individual  vortices  equals 
0.281  in  a  staggered  arrangement.   Foppl  found  that  a  pair  of 
fixed  vortices  of  suitable  strength  could  maintain  a  symmetrical 
equilibrium  position  behind  a  circular  cylinder,  provided  that 
they  are  not  subjected  to  any  asymmetric  disturbances. 

Karman' s  work  motivated  much  theoretical  work  to  extend 
his  theory  and  experimental  work  to  determine  the  validity  of 
his  stability  criterion.   The  space  ratio  has  been  generally 
confirmed  by  his  own  work  and  others  in  the  early  years .   The 
drag  coefficients  were  in  much  better  agreement  with  the  Prandtl' s 
results  than  Kirchhoff 's  even  though  it  is  still  somewhat  on  the 
low  side.   It  appeared  that  a  good  working  rule  was  available 
even  though  the  exact  mechanism  which  causes  the  vortex  shedding 
was  still  not  well  understood. 

Since  about  1950 ,  a  number  of  carefully  conducted  experi- 
ments were  designed  to  study  the  formation  of  vortex  streets .   The 
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main  effort  was  aimed  at  the  study  of  precise  dependence  of 
vortex  shedding  on  the  Reynolds  number  at  the  extreme  ranges . 
At  the  low  range  of  Reynolds  number ,  it  was  found  that  the 
"attached"  vortex  pair  postulated  by  Foppl  does  exist  but 
becomes  unstable  as  the  Reynolds  number  increases  beyond  40. 
This  wake  instability  is  responsible  for  vortex  shedding  whose 
longitudinal  spacing  differs  from  the  Kantian ' s  criterion,  and 
whose  individual  vortices  are  laminar  and  are  preserved  for 
many  diameters  downstream. 

It  was  found  that  the  shear  layers  at  the  boundaries 
of  the  wake  springing  from  both  separation  points  became  tur- 
bulent as  the  Reynolds  number  increased  to  100  and  rolled  up 
into  vortices  to  form  a  vortex  street .   The  individual  vortices 
gradually  became  turbulent  in  the  range  of  Reynolds  number  of 
150  to  300  depending  on  the  extent  of  disturbance.   Beyond 
that,  the  vortex  street  became  fully  turbulent  and  the  vor- 
tices diffused  rapidly  as  they  moved  downstream.   The  process 

remained  essentially  similar  until  Reynolds  number  reached  the 

5 
neighborhood  of  3  x  10  at  which  the  boundary  layer  itself  on 

the  cylinder  reached  the  transition  zone .   The  Strouhal  number 

remained  essentially  constant  in  this  region. 

In  the  transition  zone,  the  wake  width  became  narrow 

and  the  separation  points  moved  further  downstream  around  the 

cylinder.   This  resulted  in  a  very  rapid  increase  of  the  con- 
ventional Strouhal  number  with  Reynolds  number.   The  amount  of 
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increase  was  related  to  the  extent  of  the  motion  of  the  cylinder. 
The  increase  was  found  to  be  very  large,  if  the  cylinder  was 
held  stationary.   There  was  no  well  defined  vortex  street  in 
this  zone,  and  only  random  shedding  took  place. 

For  Reynolds  number  above  about  6  x  10  ,  it  was  observed 
that  there  was  a  sudden  decrease  in  Strouhal  number  and  a  reappear- 
ance of  a  vortex  street.   It  seems  that  proper  periodicity  does 
not  exist  in  the  turbulent  stream  above  transition.   One  may 
speak  of  a  Strouhal  number  spectrum  rather  than  a  single  value. 
The  predominant  Strouhal  number  was  determined  to  be  less  de- 
pendent on  the  Reynolds  number  and  remained  essentially  constant. 
The  magnitudes  of  the  dominant  Strouhal  numbers  were  considerably 
higher  at  the  high  Reynolds  numbers  than  the  subcritical  case. 

A  very  important  observation  made  by  Fage  and  Johansen 
was  that  the  shedding  frequency  is  related  not  to  the  diameter  of 
the  cylinder  but  to  the  width  of  the  wake.   This  led  Roshko  to 
discover  the  existence  of  a  universal  Strouhal  number  which  is 
based  on  the  wake  width,  i.e.,  the  distance  between  shear  layers 
(rather  than  the  diameter  of  the  cylinder)  and  the  wake  velocity, 
i.e.,  that  of  the  free  shear  layer  (rather  that  the  stream  velo- 
city at  infinity) .   This  number  is  essentially  constant  for  not 
only  any  Reynolds  number  but  also  for  any  common  cross-sections. 
This  constant  value  is  0.16,  see  Figure  15. 

Motivated  by  these  observations ,  Birkhof f  derived  theo- 
retically that  for  a  two  dimensional  inviscid  fluid,  the  Strouhal 
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frequency  is  inversely  proportional  to  the  geometrical  mean  of  the 
wake  width  and  the  vortex  formation  length,  i.e.,  the  average 
length  of  waving  portion  of  the  wake  before  the  shear  layers 
curl  up  to  form  a  vertex.   This  approach,  while  more  complicated, 
is  able  to  explain  many  observed  phenomena  that  contradict  the 
classical  findings. 

So  far  the  discussion  of  contemporary  investigations 
has  been  concerned  with  a  stationary  cylinder,  i.e.,  the  cylinder 
is  held  steadfast  so  that  it  cannot  move  in  response  to  the 
cross  flow  forces  resulting  from  the  vortex  shedding.   To  idealize 
the  real  situation,  in  recent  years,  experiments  have  been  con- 
ducted that  employed  spring-supported  rigid  cylinders.   It 
was  found  that,  as  long  as  the  natural  frequency  of  the  system 
is  far  removed  from  the  Strouhal  frequency,  the  cylinder  would 
vibrate  with  very  small  amplitudes.   When  they  are  near  each 
other,  linear  resonance  would  occur  to  cause  moderate  amplitude 
vibration  with  an  eddy  shedding  at  each  end  of  its  oscillation 
cycle.   In  circumstances  in  which  the  vibrating  velocity  is  com- 
parable to  the  stream  velocity,  various  nonlinear  phenomena 
occurred  more  frequently  in  cross-sections  having  smooth  edges 
and  tend  to  diminish  in  cases  involving  sharp  trailing  edges.   Also, 
the  sharp  rise  of  Strouhal  number  at  supercritical  range,  i.e., 
Reynolds  number  ranging  from  1  x  10  to  3  x  10 ,  does  not  occur 
as  long  as  the  cylinder  is  free  to  respond  to  the  cross  flow 
forces.   The  cylinder  vibrates  at  the  "correct"  Strouhal  fre- 
quencies in  such  cases. 
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When  the  cylinder  was  set  in  a  forced  vibration,  there 
was  also  a  range  of  frequencies  near  the  natural  frequency  where 
the  forced  transverse  vibrations  synchronize  or  lock-in  with 
the  vortex  frequency  to  control  the  shedding  process.   It  was 
observed  that  this  synchronized  vibration  delayed  the  initia- 
tion of  turbulence  in  the  wake  for  Reynolds  numbers  up  to  350. 
The  flow  is  essentially  two  dimensional  with  virtually  no 
appearance  of  low  frequency  irregularities  that  are  charac- 
teristics of  the  wake  behind  a  stationary  cylinder  at  low  Reynolds 
number. 

III.   MECHANISM  OF  VORTEX  SHEDDING 

The  overall  flow  pattern  around  a  circular  cylinder  with 
increasing  Reynolds  numbers  is  shown  in  Fig.  16. 

At  extremely  low  Reynolds  numbers,  say  less  than  3,  the 
flow  follows  the  Stokes  law  similar  to  a  small  particle  settling 
in  a  colloidal  solution.   The  streamlines  will  go  around  the  cy- 
linder except  the  two  which  connect  to  the  stagnation  points  and 
the  flow  pattern  is  symmetric  with  respect  to  the  cylinder  in  up- 
stream and  downstream  sides.   There  is  no  wake  in  this  range. 

At  the  Reynolds  number  of  5  to  10,  the  upper  stream  side 
remains  about  the  same  with  no  significant  change  of  pattern  but 
at  some  points  on  the  circumference  of  the  cylinder,  the  boundary 
layers  start  to  separate  in  more  or  less  symmetrical  fashion. 
The  separated  boundary  layers  beyond  the  separation  points  on  the 
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Fig.  16  Characteristic  Wake  Pattern  at  Different  Re. 
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cylinder  become  shear  layers  or  vortex  layers,  forming  a  wake 
behind  the  cylinder.   The  wake  is  a  closed  one  since  the  shear 
layers  rejoin  together  somewhere  downstream.   Inside  the  closed 
wake,  two  fixed  vortices  are  formed  immediately  behind  the 
cylinder.   These  are  the  F6ppl  vortex  pair  as  discussed  pre- 
viously.  This  pattern  is  stable  up  to  a  Reynolds  number  of 
about  40,  as  long  as  there  is  no  asymmetric  disturbance.   The 
Foppl ' s  vortex  pair  can  be  caused  to  shed  at  lower  Reynolds 
number  by  giving  the  cylinder  a  crossflow  displacement. 

The  wake  downstream  from  the  fixed  vortex  pair  becomes 
unstable  at  Reynolds  number  of  about  40  and  begins  to  shed. 
This  is  called  the  low  speed  shedding  or  shedding  due  to  wake 
instability.   The  spacing  ratio  between  the  adjacent  vortices 
is  different  from  the  value  0.281  by  virtue  of  Kantian' s  theory. 
The  wake  undergoes  a  periodic  motion  rather  than  turbulence. 

At  about  Reynolds  number  90  the  detached  shear  layers 
after  its  separation  from  the  cylinder  start   to  fold  up  into 
concentrated  vorticities .   This  is  called  the  high  speed  shedding. 
As  usual,  the  transition  zone  cannot  be  precisely  defined.   For 
a  range  of  Reynolds  numbers ,  say  between  60  to  120 ,  the  flow 
pattern  may  switch  back  and  forth  between  the  high  and  low  modes 
of  shedding  causes  irregularities  in  the  longitudinal  spacings  of 
the  vortices.   The  fluid  in  the  vortices,  however,  is  laminar  and 
the  vortex  street  persists  for  many  diameters  downstream. 

As  the  Reynolds  number  is  increased  to  the  150  to  300 
range,  the  shear  layers  beyond  the  separation  points  on  the  cylinder 
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reach  the  transition  zone  and  become  turbulent  (even  though 
the  boundary  layer  attached  to  the  cylinder  is  still  laminar) . 
The  resulting  vortices,  consisting  of  turbulent  fluid,  diffuse 
rapidly  as  they  move  downstream.   The  transition  is  complete 
at  a  Reynold  number  about  300. 

Beyond  the  transition  zone ,  the  transitional  point 
moves  closer  to  the  separation  point.   The  distance  between  them 
depends  on  the  upstream  disturbances .   The  flow  pattern  remains 
essentially  unchanged  at  least  up  to  Reynolds  numbers  of  10,000 
and  presumably  up  to  the  transition  zone  of  the  boundary  layer 
on  the  cylinder,  which  is  Reynolds  number  about  2  x  10  . 

When  the  boundary  layer  becomes  turbulent,  its  separa- 
tion point  moves  further  around  and  the  wake  width  becomes 
narrower.   If  the  cylinder  remains  stationary,  the  separation 
point  moves  back  and  forth  around  the  cylinder.   The  turbulent 
wake  is  associated  with  a  spanwise  mixed  flow.   That  is,  the 
boundary  layer  consists  of  alternate  spanwise  regions  of  laminar 
and  turbulent  flow  with  different  separation  points. 

For  Reynolds  numbers  in  the  range  of  2  x  10  and 
3  x  10  the  boundary  layer  is  undergoing  a  transition  and  the 
wake  is  disorganized.   The  formation  of  the  vortex  street  is 
suppressed.   The  energy  in  the  eddies  is  more  random  and  is 
spread  over  a  much  wider  bandwidth  of  frequency.   Beyond  Re  = 
3  x  10  ,  the  boundary  layer  is  fully  turbulent.   More  regular 
periodicity  recurs  and  a  more  or  less  regular  vortex  street  is 
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restored.   The  separation  line  on  the  rear  of  the  cylinder  becomes 

wavy,  showing  streaks  of  three  dimensional  boundary  layer  flow. 

/ 
Compared  with  the  transition  zone ,  the  excitation  in  full  turbulent 

flow  is  more  regular  and  it  may  be  called  "narrow-band  random" 
or  quasi-periodic.   At  Reynolds  number  of  10  x  10  ,  the  separation 
line  becomes  straightened  and  moves  slightly  upstream.   Mean- 
while the  disturbance  at  the  rear  face  of  the  cylinder  gradually 
disappear  at  higher  Reynolds  numbers.   This  suggests  that  at  even  higher 
Reynolds  number,  the  quasi-periodic  pattern  should  prevail. 
There  is  no  evidence  to  verify  this,  partially  because  of  the 
extreme  difficulties  in  conducting  experiments  in  a  high  Reynolds 
number  regime  such  as  the  larger  throat  area  and  high  pressure 
requirements . 

The  boundary  layers  after  separating  from  the  cylinder 
become  free  shear  layers .  After  several  diameters  the  shear 
layers  roll  up  to  become  the  vortex  street.   This  stretch  is 
called  the  formation  length.   The  flow  pattern  is  very  sensi- 
tive to  the  disturbance  there.   The  flow  pattern  will  be  changed 
if  fluid  is  injected  into  the  wake.  >  This  can  be  easily  done 
through  the  body  itself.   This  technique  has  been  given  the  term 
base  bleeding.   For  low  bleeding  rates,  the  wake  narrows  down 
as  if  the  body  were  elongated  into  the  wake.   At  higher  bleeding 
rates  the  wake  is  widened.   For  very  thin  cylinders,  say  h   mm 
in  diameter,  the  pattern  is  disturbed  even  by  insertion  of  a  hot 
wire  in  the  shear  layer. 
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There  is  evidence  which  indicates  that  the  free  shear 
layers  from  either  side  of  the  wake  interfere  with  each  other  in 
the  formation  length.   If  a  thin  partition  is  placed  along  the 
centerline  of  the  wake,  downstream  of  the  cylinder,  the  pattern 
of  flow  will  be  changed.   This  technique  has  been  given  the  term 
"splitter  plate".   If  the  length  of  a  splitter  plate  is  long, 
the  vortex  formation  is  totally  prevented.   For  a  short  pldte, 
placed  near  the  cylinder,  the  effect  is  similar  to  low  base  bleed- 
ing and  the  vortices  form  downstream  from  the  splitter  plate. 
This  downstream  shift  of  formation  location  tends  to  narrow  the 
wake.   On  the  other  hand,  when  the  plate  is  placed  far  downstream 
beyond  the  formation  length,  the  vortices  form  in  normal  way  and 
no  change  in  flow  pattern  is  observed.   Pressure  measurement 
showed  that  the  minimum  pressure  was  moderated  by  the  presence  of 
the  splitter  plate  and  its  location  moved  further  downstream. 

The  flow  pattern  is  also  sensitive  to  the  surface  condition 
of  the  cylinder  such  as  roughness.   The  cause  may  be  attributed 
to  interference  with  the  separation  line .   Experiment  showed  that 
by  placing  trip  wires  on  the  surface  of  the  cylinder,  the  separa- 
tion line  can  be  significantly  altered.   It  was  found  that,  re- 
gardless of  the  disturbance,  a  vortex  street  can  exist  as  long  as 
a  more  or  less  straight  separation  line  can  be  maintained  on  the 
cylinder  surface.   Vortex-shedding  can  be  virtually  inhibited  by 
a  crisscross  arrangement  of  trip  wires,  for  which  a  straight  separa- 
tion line  is  impossible. 
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The  width  of  vortex  street  is  not  exactly  constant. 
Careful  observations  show  that  the  transverse  spacing  between 
the  vortex  rows  first  increases  downstream  from  the  cylinder, 
then  attains  a  maximum  and  finally  decreases  again  before  the 
street  becomes  unstable  and  disintegrates. 

The  growing  vortex  is  fed  by  circulation  from  the  shear 
layer  until  strong  enough  to  begin  rolling  up  and  drawing  the 
opposite  shear  flow  across  the  wake.   The  approach  of  vorticity 
of  opposing  sign  cuts  off  further  circulation  to  the  growing 
vortex,  which  is  then  shed  and  moves  downstream  . 

If  the  section  has  sharp  corners,  the  flow  separation 
will  occur  there,  for  all  different  Reynolds  numbers.   Thus,  the 
flow  patterns  will  be  less  dependent  on  the  Reynolds  numbers. 
However,  the  section  having  sharp  corners  will  have  stronger 
vortex  excitation  and  will  be  more  unstable  than  a  smooth  sec- 
tion.  This  is  explained  by  the  better  correlation  of  disturbance 
along  the  span. 

Another  important  sectional  shape  parameter  is  the 
streamwise  length  and  the  shape  of  the  "afterbody"  -  which  is 
defined  as  the  part  of  the  section  downstream  of  the  separation 
points,  or  the  sharp  edge.   For  example,  the  D-section  with  the 
flat  face  upstream  will  experience  large  excitation  while  with 
the  flat  face  downstream  will  experience  a  negligible  excitation. 
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IV.   STROUHAL  NUMBERS 

It  is  well  known  that  the  Strouhal  Number  is  shape  dependent.   The 
most  widely  reported  Strouhal  Number  is  associated  with  a  circular 
cylinder.   The  Strouhal  Number  for  rectangular  sections  are  also  fairly 
well  studied.   However,  for  other  sections,  especially  structural 
shapes,  there  is  only  scant  information  available  in  literature. 

The  classical  Strouhal  Number  for  a  cylinder  placed  normal  to 
the  current  is  defined  as: 

st  =  M 

V 

where  St  =  Strouhal  Number 

f  =  frequency  in  Hertz 

d  =  diameter  of  cylinder 

•V  =  speed  of  stream  current 
1.   Circular  Cylinders 

For  a  smooth  circular  cylinder  the  relation  between  the  Strouhal 
Number  and  the  Reynolds  number  is  shown  in  Fig.  17   .   It  is  seen  that 
it  is  roughly  divided  into  six  regions.   Strouhal ' s  phenomenon,  as 
explained  by  Rayleigh,  that  Strouhal  Numbers  for  a  thin  wire  decreases 
slightly  with  a  rising  temperature,  (hence,  with  an  increase  in 
Reynolds  number)  is  true  only  in  the  region  of  Turbulent  Vortex 
Street-Laminar  Boundary  Layer.   For  low  Reynolds  numbers,  which  com- 
prise the  Transition  range,  Laminar  Vortex  Street  range,  and  Wake 
Instability  Vortex  Street  range,  the  Strouhal  Number  decreases  rapidly 
with  the  decrease  of  Reynolds  numbers.   By  an  empirical  fitting,  it 
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was  found  that  the  following  expressions  hold  true  approximately: 

St  =  0.195  -  3.92/Re  (Rayleigh,  Tyler) 
but  a  better  fit  is  given  by 

St  =  0.212  -  4.5/Re  (Roshko) 
for  Reynolds  number  up  to  about  250, 
and  St  =  0.212  -  2 . 7/Re  (Roshko) 

'for  Reynolds  numbers  in  -the  range  of  250  to  2000.  •   The  detailed 
relationship  is  shown  in  Fig.  18  .   Beyond  Re  =  2000,  the  flow  is  the 
Turbulent  Vortex  Street  range,  Strouhal  Numbers  appear  to  decrease 
from  0.21  gradually  to  about  0.18,  when  the  flow  reaches  critical  con- 
dition, and  enters  the  incipient  transition  flow  .   This  range  was 
studied  by  Strouhal. 

The  exact  Reynolds  number  at  which  the  critical   flow   starts 
cannot  be  explicitly  stated.   It  can  be  as  high  as  4  x  10  ,  depending 
on  the  disturbances.   The  disturbances  can  originate  from  the  tur- 
bulent flow  upstream  or  from  the  irregularities  on  the  surface  of  the 
cylinder.   Extreme  roughness  on  the  surface  of  the  cylinder  will 
greatly  advance  the  incipient  transition  flow. 

In  the  transition  range,  the  Strouhal  Numbers  can  assume  a  wide 
range  of  values ,  depending  principally  on  the  support  condition  of 
the  cylinder,   (See  Fig.  19  ),  if  the  cylinder  is  held  stationary 
(i.e.,  motionless) ,  the  Strouhal  Number  can  be  very  high,  with  a  reported 
value  as  much  as  0.5.   If  the  cylinder  is  free  to  respond  to  exci- 
tations, a  Strouhal  Number  may  be  as  low  as  0.2,  which  was  observed 
in  cases  of  large  smoke  stacks.   Various  intermediate  values  have 
also  been  recorded  and  some  observers  insist  that  the  shedding  is 
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aperiodic.   Observations  of  spectra  of  velocity  fluctuations  indicated 
that  the  shedding  extended  over  a  broad  band  of  frequencies.   Each 
frequency  in  the  band  carries  a  portion  of  the  energy  from  the  wind 
stream.   The  frequency  which  carries  the  largest  portion  of  the  energy 
is  the  dominant  frequency  at  which  the  cylinder  presumbly  vibrates. 

The  fact  that  the  energy  is  distributed  over  a  broad  band  of 
frequencies  in  a  transition  range  is  of  particular  interest  in  the 
structural  design  criterion  concerning  wind-induced  vibrations.   It 
is  clear  that  the  energy  spectrum  would  be  flatter  over  a  broad 
band  of  frequencies  and  hence  the  root  mean  square  energy  for  any 
given  frequency  would  be  smaller.   It  is  for  this  reason  that,  in 
a  good  design  practice,  a  suitable  adjustment  of  design  parameters 
should  be  made  to  keep  Reynolds  numbers  in  the  transition  range. 

The  shape  of  energy  spectrum  depends  on  the  shape  of  the  speci- 
men and  upstream  turbulence,  etc.,  see  Fig.  20.     It  is  most  probably 
dependent  also  on  the  support  condition  of  the  cylinder.   The  spectra 
would  shift  substantially  downward  in  the  frequency  scale  if  the 
cylinder  is  permitted  to  vibrate  with  the  excitations.   In  the  current 
literature,  investigations  are  limited  to  a  small  range  of  Reynolds 
numbers,  and  are  applicable  only  to  stationary  circular  cylinders. 

In  the  Trans-critical  range ,  pseudo-regular  or  quasi-periodic 
shedding  resumes  and  the  frequency  response  is  narrow-banded  (but  not 
as  narrow  as  in  the  subcritical  range) .   The  Strouhal  Number  is  grad- 
ually on  the  rise  with  the  increase  of  Reynolds  number.   This  rise 
is  due  to  the  downstream  shift  of  the  separation  points  of  the  shear 
layer  on  the  cylinder  wall.  A   representative  figure  for  the  Strouhal 
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Number  is  0.27  for  Reynolds  number  =  5.0  x  10  ,  and  slightly  higher  for 
higher  Reynolds  numbers. 
2.   Rectangular  Section 

Two  factors  are  important  for  rectangular  sections :  sharpness  of 
the  corners  and  the  length  of  the  after-body.   For  a  truly  sharp 
corner,  the  separation  will  always  occur  on  the  corner  point,  indepen- 
dent of  the  Reynolds  number.   Strouhal  Numbers  in  this  case  will  be 
nearly  independent  of  Reynolds  numbers.   For  a  rectangular  section 
which  has  round  corners ,  there  will  be  some  downstream  receding  of 
the  separation  points  at  increasing  Reynolds  numbers ,  with  a  corres- 
ponding increase  of  Strouhal  Numbers  as  a  consequence  of  it,  although 
this  increase  is  more  moderate  in  comparison  to  the  case  for  a  circular 
cylinder. 

The  other  important  sectional  parameter  is  the  streamwise  length 
and  shape  of  the  part  of  the  section  downstream  of  the  separation 
points.   This  so-called  after-body  is  the  part  of  the  structure  most 
directly  exposed,  to  the  vortex  excitation.   A  long  after-body  will 
interfere  with  the  wake  vortex  formation,  much  like  a  splitter  plate 
and  can  also  cause  the  separated  flow  to  reattach  to  its  surface  at  a 
point  further  downstream.   The  flow  may  separate  again  after  reattach- 
ment for  a  very  long  afterbody,  and  the  vortex  formation  takes  place 
in  the  wake  of  the  second  separation.   In  either  case,  the  Strouhal 
Number  is  decisively  influenced  by  this  parameter.   Fig.  21   shows  the 
variation  of  Strouhal  Numbers  with  the  after-body  length  ratio.   It  is 
seen  that  the  Strouhal  Number  assumes  a  constant  value  of  about  0.14 
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except  when  h/d  is  in  the  range  of  1.5  to  3  where  Strouhal  Number  can 
be  as  low  as  0.06. 
3.   Other  Shapes 

As  the  radius  of  the  round  corner  in  the  rectangular  section 
increases  to  a  limit,  the  section  becomes  an  ovaloid,  which  resembles 
an  ellipse.   Ovaloidal  as  well  as  elliptical  sections  have  been  in- 
vestigated by  a  few  researchers  for  a  limited  range  of  Reynolds  num- 
bers.  The  results  are  shown  in  Fig.  22.    It  is  seen  that  the  results 
are  in  general  agreement  with  what  might  have  been  expected  from  the 
extrapolation  from  the  results  for  rectangular  sections.   The  support- 
ing condition  is  not  given  but  the  rather  high  Strouhal  Numbers  in 
the  transition  range,   seem  to   indicate  that  the  specimens  were 
probably  held  stationary  during  the  tests.   If  the  specimens  were 
allowed  to  vibrate  with  the  excitations ,  the  dominating  Strouhal 
Numbers  should  not  differ  much  from  those  for  subcritical  in  the  trans- 
critical  range.   Further  experiments  are  needed  to  verify  this  con- 
jecture. 

Scattered  results  for  other  miscellaneous  shapes  of  interest  to 
aviation  and  space  technology  such  as  diamond,  triangle,  and  crescent 
are  available  in  the  literature.   The  results  will  not  be  presented 
here,  since  they  are  rather  incomplete  and  the  sections  themselves 
have  limited  applications  to  structural  design.   It  is  sufficient  to 
say  the  Strouhal  Number  is  generally  in  the  range  of  0.15  to  0.2  and 
increases  gradually  with  the  increase  of  Reynolds  number  up  to  the 
transition  range.   If  the  member  is  held  stationary,  the  Strouhal 
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Number  will  exhibit  a  characteristic  jump  to  around  0.4  to  0.6,  similar 
to  the  case  of  cylinders  and  rectangular  sections,  where  the  Reynolds 
numbers  fall  in  the  transition  zone.   Otherwise,  the  cylinder  will 
vibrate  with  excitations  at  the  "correct"  frequencies. 

Of  particular  interest  to  the  bridge  designers  are  the  structural 
rolled  shapes,  which  includes  HP,  I,  Channels,  Angles,  Wide-Flanges, 
and  Built-Up-Sections .   The  best  source  of  information  on  Strouhal 
Numbers  for  rolled  shape  is  ASCE  Task  Force  on  Wind  Force  Committee 
Report,  which  contains  a  comprehensive  table  (Fig. 23   ) .   It  is  seen 
that  a  variety  of  structural  shapes  and  built-up  sections  are  included, 
However,  with  the  exception  of  the  cylinder,  no  indication  of  Reynolds 
numbers  is  given.   From  the  discussion  in  the  previous  section,  all 
the  structural  shapes  are  sharp- cornered.   Since  the  separation  points 
of  the  shear  layers  must  be  fixed  on  the  sharp  corners ,  the  depen- 
dence of  Strouhal  Numbers  on  the  Reynolds  numbers  is  probably  not  very 
strong. 

The  variation  of  Strouhal  Number  between  different  cross-sectional 
shapes  (regardless  the  wind  direction)  is  not  very  large.   For  example, 
among  all  the  sections  which  possess  two  axes  of  symmetry,  the  Strouhal 
Number  is  a  minimum  for  HP  sections  with  a  value  of  0.12  and  is  maxi- 
mum for  twin  channel  sections  with  a  value  of  0 . 15 ,   For  unsymmetric 
sections,  the  Strouhal  Number  can  vary  between  0.114  to  0.168.   The 
only  other  measurement  known  to  us  is  by  Narita  for  H  shapes  with 
deptsh/width  ratios  of  1/2  and  3/4.   The  reported  Strouhal  Numbers 
were  0.115  and  0.12,  respectively.   These  are  in  reasonably  agreement 
with  data  given  in  Fig.  23. 
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From  the  above  discussion,  it  appears  to  be  reasonable  to  esti- 
mate the  Strouhal  Number  of  HP  sections  to  be  0.12  for  engineering  pur- 
poses.  In  view  of  the  fact  that  the  most  useful  range  of  Reynolds 
numbers  should  be  in  the  transition  range  in  which  the  excitation  fre- 
quencies are  broad-banded,  a  crude  estimate  of  the  Strouhal  number  should 
not  lead  to  serious  consequences. 

On  the  other  hand,  it  is  regrettable,  or  even  embarrassing,  that 
members  of  great  engineering  application  such  as  HP  sections  had  not 
stimulated  more  than  cursory  attempts  in  the  determination  of  their 
Strouhal  frequencies.   A  more  comprehensive  program  to  determine  the 
Strouhal  Numbers  of  commonly  used  structural  shapes  should  be  a  worth- 
while undertaking. 
4.   Effect  of  Wind  at  an  Obligue  Direction 

It  is  of  great  practical  interests  to  determine  the  Strouhal 
Number  when  the  wind  direction  is  at  an  angle  other  than  90  degrees 
with  the  longitudinal  axis  of  the  circular  cylinder.   This  type  of 
investigation  has  become  popular  during  the  last  decade.   The  range 
of  yaw  angles  covered  was  very  extensive  but  the  range  of  Reynolds 
numbers  is  limited  to  below  2  x  10  . 

It  is  expected  that  the  flow  would  have  both  longitudinal  and 
transverse  components.   If  the  shear  layer  separation  is  unaffected 
by  the  presence  of  spanwise  flow  and  the  vortex  street  is  undisturb- 
ed by  the  yaw  angle,  vortices  will  shed  at  a  frequency  corresponding 
to  the  normal  velocity  component,  i.e., 

St   =  St  Sina 

a 
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where  a  =  angle  between  the  flow  and  the  longitudinal  axis  of  the  sta- 
tionary cylinder.   Observations  seem  to  confirm  this  assumption.   Simi- 
lar results  were  observed  for  yawed  music  wire,  which  can  freely  res- 
pond to  the  excitation.   When  a  is  less  than  18  degrees,  there  is  some 
scatter  of  frequencies  with  more  broad  band  width.   Fig.  24   shows- 
the  comparison  between  the  test  results  and  the  computed  frequencies 
based  on  the  Sine  Law.   It  is  seen  that  the  agreement  is  rather  good, 
at  least  for  low  Reynolds  numbers .   Very  little  is  known  about  the 
yaw  effect  at  high  Reynolds  numbers.   There  is  no  reason  to  expect 
that  the  results  would  be  different  from  that  for  lower  Reynolds 
numbers . 

When  the  wind  is  applied  at  an  angle  of  inclination  to  the 
width  of  the  body,  the  effect  depends  on  the  length  of  the  after-body 
as  discussed  before.   The  lift  coefficient  would  be  very  sensitive 
to  the  angle  of  inclination.   Also,  the  flow  pattern  may  be  very 
different  and  sometimes  the  galloping  vibration  is  liable  to  occur 
at  relatively  high  angle  of  inclination.   This  discussion,  however, 
is  beyond  the  scope  of  this  report. 
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V.   SUPPRESSION  OF  VORTEX  SHEDDING 

There  are  three  basic  methods  to  mitigate  the  vortex-induced 
vibrations : 

(1)  Aerodynamic 

(2)  Mechanical 

(3)  Analytical 

All  these  methods  have  been  used  to  varying  degrees  of  success. 
Systematic  procedure  or  devices  have  been  derived  from  these  basic 
methods  and  used  in  practical  cases.   The  results  were  satisfactory  or 
at  least  encouraging.   The  following  sections  will  treat  each  method 
in  detail. 
1  ,  Aerodynamic  Method 

The  aerodynamic  method  of  suppressing  the  vortex  induced  vibra- 
tion entails  the  disruption  of  the  vortex  forming  mechanism.   From  p. 46 
Qiapter  III  it  is  stated  that  alternating  shedding  of  vortices  can 
occur  only  in  a  bluff  body  in  which  the  boundary  layers  can  separate 
to  free  shear  layers.   To  have  a  strong  vortex  street,  the  separation 
points  must  be  lined  up,  preferably  in  a  straight  line  perpendicular 
to  the  current  direction,  so  that  a  spanwise  coherence  is  effected. 
Furthermore ,  the  shear  layer  must  have  an  adequate  vortex  formation 
length  to  roll  up.   The  area  behind  the  bluff  body  must  be  maintained 
as  a  fluctuating  but  essentially  low  pressure  region.   If  any  or  all 
of  these  conditions  are  disrupted,  the  vortex-shedding  vibration  will 
be  greatly  attenuated  or  even  stopped  altogether. 
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The  following  devices  were  designed  for  the  purpose  of  disrupt- 
ing some  of  these  conditions : 

A.  Streamlined  Fairings,  which  change  the  bluff  bodies  into  stream- 
lined bodies. 

B.  Roshko's  Splitter  Plate  or  Baird's  Saw-Toothed  Paneling,  which 
interferes  with  the  fluctuation  of  the  localized  low  pressure  re- 
gion in  the  aft  of  the  cylinder. 

C.  Perforations  in  the  Bluff  Body,  which  eliminate  the  low  pressure 
region  in  the  aft  by  means  of  base  bleeding. 

D.  Perforated  shroud  over  the  Bluff  Body,  which  interferes  with  both 
the  separation  of  shear  layer  and  the  low  pressure  region. 

E.  Trailing  Ribbons  which  interfere  with  the  shear  layers  in  the  vor- 
tex formation  length. 

F.  Helical  Strakes  or  Spiral  Wrapping  which  disrupt  the  spanwise 
coherence  of  vortex  shedding. 

The  design  of  these  devices  is  at  present  beyond  the  capability 
of  analytical  means  and  must  be  done  by  trial  and  error  method  in  a 
wind  tunnel.   The  following  paragraphs  summarized  the  main  results  and 
design  criteria. 

(A)   Streamlined  Fairing 

The  streamlined  fairings  have  been  in   widespread  usage 
for  many  years  and  considerable  research  effort  has  been  directed  to 
assess  their  hydrodynamic  performance.  In  general,  a  properly  designed 
fairing  can  be  effective  as  long  as  the  wind  direction  remains  unchanged. 
Also,  the  long  ' afterbody'  produced  by  a  well-proportioned  fairing  is 
amenable  to  galloping. 
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(B)   Splitter  Plate 

A  splitter  plate  is  a  thin  partition  placed  along  the  cen- 
ter line  of  the  wake  downstream  of  a  cylinder.,   It  appears  that  it  caus- 
es the  two  shear  layers ,  in  the  formation  length ,  to  be  isolated  from 
each  other  and  broken  down  independently.   It  also  has  the  effect  of 
increasing  the  back  pressure  and  thus  reducing  the  pressure  gradient 
from  the  separation  points  to  the  wake.   The  indirect  result  is  a 
reduction  of  drag  coefficient  of  the  cylinder.   In  the  subcritical 
region  of  Reynolds  numbers ,  the  splitter  plate  fcrces  the  vortices 
to  form  further  downstream  along  its  length.   This  downstream  shift 
of  formation  position  is  accompanied  by  a  decrease  of  the  Strouhal 
frequency.   The  effective  length  of  the  splitter  plate  was  found  to  be 
about  4  to  5  cylinder  diameters.   Limited  tests  on  splitter  plate  for 
supercritical  range  indicated  that  the  effectiveness  was  much  less 
noticeable.   The  splitter  plates  are  most  effective  near  antinodes  of 
a  vibrating  flexible  cable.   For  economic  reasons,  the  splitter 
plates  need  not  be  continuous  along  the  entire  length  of  the  cylinder, 
but  rather,  they  can  be  placed  intermittently  along  the  length,  each 
segment  being  about  10  to  20  diameters  long.   Also,  Baird  found  that 
segment  of  plates  of  varying  widths,  say  triangular  shaped,  (a  bank 
of  the  plates  would  have  an  appearance  of  saw  teeth)  were  more  effective 
than  those  of  constant  width.   Alexander  found  that  random  arrangement 
of  longitudinal  fins  was  also  an  effective  design. 
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(C)  Perforations 

Perforations  (Fig.  25)  are  holes  of  suitable  sizes  and 
shapes  on  the  bluff  body,  which  allow  a  sufficient  amount  of  addi- 
tional fluid  directly  into  the  localized  low  pressure  region  in  the 
aft.   It  was  found  that  porosity  of  0.5  was  effective  in  inhibiting 
the  vortex  shedding.   The  size,  shape  and  distribution  of  the  holes  are 
apparently  unimportant  in  the  performance  of  porosity,  only  the  total 
hole  area  relative  to  the  total  frontal  area  matters.   For  economical 
reasons ,  it  is  easier  to  provide  a  few  large  holes  than  a  large  num- 
ber of  small  holes.   Also,  the  curvature  of  the  hole  boundary  should 
be  as  flat  as  possible  to  minimize  the  effect  of  stress  concentration. 
It  appears  that  the  ovaloidal  holes  equally  spaced  on  the  surface 
through  the  body  would  be  a  good  choice.   This  arrangement  has  been 
used  very  effectively  in  the  box  members  in  trussed  bridges  for  weight 
reduction  purposes .   The  same  design  can  be  used  also  for  vibration 
suppression. 

(D)  Shroud 

Shroud  is  a  perforated  cylindrical  shell  enclosing  the 
bluff  body  with  a  small  diameter  gap, (Fig.  26).   It  is  found  that  the 
amount  of  perforation  relative  to  the  total  area  should  be  about  20 
to  30  percent  and  the  gap  should  be  about  0.12  times  the  diameter. 
This  method  has  been  successfully  applied  to  anti-corrosion  anodes  of 
cooling  water  ducts ,  multi-flue  chimney  stacks ,  and  square-section 
towers . 
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(E)  Trailing  Ribbon 

Trailing  ribbons  are  thin  plastic  film  strips  attached  to 
the  bluff  body  so  that  they  flutter  in  the  vortex  forming  length 
downstream  from  the  bluff  body.   The  fluttering  motions  of  the  ribbons 
interfere  with  the  rolling  up  of  the  shear  layers  to  prevent  them  from 
forming  a  vortex  street.   It  has  been  used  successfully  in  marine  appli- 
cations, to  suppress  the  strumming  of  a  towed  array.   It  was  found 
there  that  well-designed  ribbons  reduced  both  the  frequency  and  the 
amplitude  of  vibration,  at  least  in  water, for  Reynolds  numbers  in 
the  subcritical  range.   It  is  presumed  that  similar  effect  could  also  be 
realized  for  a  bluff  body  in  wind.   The  optimum  ribbon  configuration 
was  found  to  be  6  and  10  d  long*,  by  2  d  wide  spaced  1  and  2  d  apart. 
In  an  oblique  stream,  the  ribbons  have  a  tendency  to  lay  flat   along 
the  longitudinal  axis  of  the  cylinder  if  any  appreciable  spanwise 
stream  flow  is  present.   In  such  circumstances  the  effectiveness  of 
the  trailing  ribbons  would  be  greatly  impaired. 

(F)  Helical  Strakes 

Helical  strakes  (Fig.  27)  are  ridges  wrapped  around  a 
cylindrical  or  nearly  cylindrical  body  to  trip  the  separation  points 
on  the  cylinder.   Due  to  the  helical  pattern  of  the  strakes,  the 
separation  points  would  be  expected  to  vary  along  the  axis  of  the 
cylinder.   The  vortex  shedding  would  become  incoherent  spanwise  and 

*  "d"  denotes  the  diameter  of  the  cable 
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the  fluctuating  lift  force  thus  would  be  greatly  weakened.   Many  factors 
affect  the  performance  of  the  strakes ,  which  are:  height,  width,  and 
shape  of  the  strakes,  and  the  pitch  of  the  helical  winding.    Tests  by 
Scruton,  Woodgate ,  Fabula  and  others  show  no  conclusive  results.   It 
is  generally  agreed,  however,  that  the  ridge  height  should  be  0.1  d 
or  more  to  be  effective. 
2  .  Mechanical  Method 

There  are  two  main  types  of  devices  to  minimize  vibrations  me- 
chanically:  Tie-downs  and  Damping.   Tie-downs  are  struts  which  con- 
nect the   expected  anti-nodes  of  the  member  to  some  rigid  supports. 
The  amount  of  force  taken  by  the  strut  is  dependent  on  the  aerody- 
namic lift  force   of  the  member,  which  is  usually  rather  small  in 
terms  of  structural  strength. 

Dampers   are  a  mass -spring  combination  which  has  same  natural  fre- 
quencies as  one  of  the  natural  frequences  of  the  structure  itself. 
The  energy  is  transferred  from  the  structure  to  the  damper,  and  is 
dissipated  by  friction  with  viscoelastic  pads,  synthetic  rubber 
stemmed  damper,  etc.   Spring  can  be  leaf,  helical,  or  torsion  types, 
and  the  position  of  the  weight  is  adjustable.   The  damper  usually 
needs  frequent  maintenance  and  is  not  very  effective.   A  Stockbridge 
damper  (Fig.  28)  consists  of  a  double  cantilever  attached  to  a 
vibrating  structure .   The  cantilever  is  tuned  to  vibrate  with  the 
same  frequency  as  the  vortex  shedding  frequency  but  is  ir-radians  out 
of  phase.   The  inertial  force  is  therefore  cancelled  out  this  way.  It 
is  not  a  true  damper  in  the  ordinary  sense.   It  is  generally  used  for 
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light  structures  such  as  a  wire  but  the  same  principle  can  be  used  for 

very  heavy  structures , even  to  restrain  a  ship  from  excessive 

rolling. 

3.   Analytical  Method 

It  is  recommended  by  some  that  the  violent  vibration  due  to 
resonance  can  be  designed  out  by  setting  the  natural  frequency 
corresponding  to  a  wind  velocity  much  above  the  prevailing  wind 
velocity  near  the  site  of  the  structure.   This  is  not  always 
possible  and  is  very  wasteful  even  if  it  is  possible. 

A  more  practical  idea  is  to  design  the  member  such  that  the 
Reynolds  numbers  corresponding  to  the  prevailing  wind  speeds  be 
kept  within  the  supercritical  regime  (random  shedding  range, 
Fig.  17) .   We  recall  that  the  excitation  frequency  within  the 
transition  zone  is  broad-banded,  (Figs ,19  &  20) .   It  is  clear  then 
that  the  dominating  frequency  would  carry  a  relatively  small 
share  of  total  energy  of  the  wind.   In  such  cases,  the  excursion 
and  bending  moment  are  expected  to  be  relatively  small  (see 
example  5 ,  Part  C,  Chapter  III) .   The  analytical  method  can  be  a 
very  useful  tool  when  it  is  used  in  conjunction  with  some  or  all 
of  the  other  methods . 
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VI.    DESIGN  DATA 

(1)   End  -  Fixity  Coefficients 

In  the  determination  of  the  natural  frequencies  and  mode  shapes 
of  a  bridge  member,  it  is  necessary  to  know  the  degree  of  fixity  of 
that  member  at  its  ends.   This  can  be  accomplished  easily  as  follows: 
The  member  should  be  considered  supported  at  each  end  by  springs  which 
represent  the  action  of  the  surrounding  members  (See  Fig.  below.)   If 
the  spring  constant  k  can  be  determined  for  each  spring,  then  the  degree 
of  fixity  can  be  determined.   From  the  standard  theory  of  structures 

it  is  known  that,  for  a  unit  moment 


1  imposed  on  an  end  of  the  member ,  the 

angle  of  rotation  that  results  in  the 
member  at  that  end  is  given  by: 

n 


I. 

Member  A  4  E   ^  — 

1  H 

I. 

where  £  — —  equals  the  summation  of  the  stiffness  of  adjacent 

i   L. 

l 

members  of  the  particular  joint  containing  the  end  of  the  member  being 
analyzed  n  is  a  constant  which  depends  on  the  continuity  of  the  member  in 
consideration  relative  to  the  joint.   If  the  member  is  continuous  through 
the  joint  (e.g.  member  A)  then  n  =  2;  this  is  commonly  the  case  for  a 
lower  chord  in  a  truss.   If  the  member  is  terminated  at  the  joint  (e.g., 
member  B)  then  n  =  1;   this  is  commonly  the  case  for  a  hanger  of  a  bridge 
deck. 

The  spring  constant  at  an  end  of  a  member  is  then 

k-i.-*L    z^i 

6    n   i   l± 
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and  the  nondimensional  spring  constant,  namely,  the  "end  fixity 

condition"  is 

£  _i 

kL    4   i  L. 


a  = 


El 


I 
L 


where  L  =  Length  of  the  member  in  consideration, 

I  =  Moment  of  inertia  of  the  member  in  consideration. 

Specifically,  from  the  definition  of  n  stated  above,  the  end-fixity 

conclusions  are 

4  I    Ii/L. 
a  =   for  hangers  (member  B) 


a 


I/L 


2  Z  1±/L± 

l 


for  lower  chord  (member  A) 


I/L 


The  above  expressions  show  that  the  end-fixity  condition  a 
has  the  following  physical  interpretation:  it  is  equal  to  four  times 
the  ratio  of  the  joint  stiffness  and  the  flexural  stiffness  of  the 
member  in  question. 


Member  B 
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(2)   Damping  Coefficients 

The  damping  in  a  bridge  member  is  due  to  a  combination  of  struc- 
tural (material),  viscous  and  aerodynamic  effects.   The  damping  co- 
efficient is  usually  expressed  in  the  logarithmic  decrement,  which 
varies  with  amplitude,  mode  (flexural  or  torsion)  and  material.   For 
members  which  have  considerable  flexural  strength,  such  as  a  struc- 
tural shape  or  a  built-up  section,  the  aerodynamic  damping  is  negli- 
gible and  the  damping  is  primarily  due  to  structural  damping  at  low 
excursions  and  due  to  viscous  damping  at  high  amplitudes. 


Whitebread  suggested  values  of  logarithmic  decrement  for  various 
types  of  construction : 

Type  of  Construction  Logarithmic  Decrement 

Lightweight  aluminum  0.005  to  0.01 

Welded  steel  (unlined)  0.01  to  0.05 

Welded  steel  (lined)  0.03  to  0.07 

Concrete  (in  situ)  0.05  to  0.10 

Concrete  (precast)  0.07  to  0.15 

Narita  reported  some  experimental  damping  factors  for  three 
steel  bridges  in  Japan,  and  one  in  the  U.S.S.R.   The  logarithmic  decre- 
ment values  for  the  hangers  range  from  0.008  to  0.032  for  built-up 
plates  and  from  0.07  to  0.014  for  pipes. 
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Selberg  reported  that  the  logarithmic  decrement  for  steel 
bridges  ranges  from  0.02  to  0.05. 

Sachs  reported  that  the  logarithmic  decrement  for  cables  is 
about  0.002  which  is  much  lower  than  built-up  section  and  the  pipe. 
Table  2  shows  damping  coefficients  of  chimneys, towers,  etc.,  appeared 
in  literature. 
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(3)   Lift  Coefficients 

Since  Lord  Rayleigh  observed  that  the  vortex-shedding  vibrations 
were  affected  in  a  plane  perpendicular  to  the  stream  flow,  it  follows 
that  there  must  exist  a  fluctuating  lift  force  on  the  specimen. 
Early  attempts  to  measure  this  lift  force  were  made  by  Schwabe, 
Bingham  and  others.  Many  of  these  early  experiments  were  handicapped 
by  inadequate  instrumentation  and  poor  test  procedure.  Their  results 
are  widely  scattered  and  must  be  treated  with  caution. 

The  first  controlled,  systematic  measurement  of  lift  acting  on 
a  cylinder  at  high  Reynolds  numbers  were  carried  out  by  Fung  in  1960. 
The  tests  were  duplicated  by  Schmidt  four  years  later.   The  rms  lift 
coefficients  from  these  sources  are  shown  in  Fig.  29.   It  is  seen 
that  one  is  nearly  three  times  as  high  as  the  other.   Schmidt 
attributed  the  difference  to  the  fact  that  Fung's  model  had  minor 
surface  irregularities  while  his  was  highly  polished.   Keefe  also 
found  that  gap  ,  vents  and  tip  shapes  at  the  ends  of  the  cylinder  had 
a  tendency  to  increase  the  unsteady  lift  force.   Even  a  trip  wire 
or  a  thin  tape  bonded  longitudinally  on  the  cylinder  would  do  the 
same.   It  appears  that  the  lift  force  is  affected  by  the  spanwise 
coherence  of  the  shedding  mechanism.   In  fact,  the  lift  coefficient 
is  also  closely  related  to  the  drag,  pressure  and  Strouhal  number. 
Since  the  Strouhal  number  in  post-critical  regime  is  strongly  affected 
by  the  movement  of  the  cylinder,  the  lift  coefficient  must  be  also. 
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The  relation  between  the  rms  and  peak  lift  coefficients  are 
shown  in  Fig.  30,  as  reported  by  Fung.  It  is  seen  that  the  peak 
lift  is  about  twice  of  the  rms  value. 

In  the  limited  literature  available  the  lift  coefficient  CL 
was  reported  by  Narita  who  estimated  it  to  be  approximately  1.15 
for  3  different  types  of  H  shapes. 
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PART  B 
MATHEMATICAL  THEORY 
I.   NATURAL  FREQUENCY  AND  NORMAL  MODES 


(1)   Flexural  Vibrations 

The  small  amplitude ,  free ,  transverse  vibration  of  a 
beam  is  governed  by 

d   (EI   d_Y  )  -  (^  /y  =  0 


dX2"      dX7 
where  Y  =  transverse  displacement 

EI  =  flexural  rigidity 

PT  =  mass  density 

A  =  cross-sectional  area 

a)  =  circular  frequency. 

If  axial  tension  T  is  imposed,  the  equation  becomes 

_4  (EI^J  )  -  T  ^     -pA^Y  =  °  (1) 

dx2^    dsr  ax1 

For  prismatic  beams,  EI  =  constant.   Nondimention- 
alizing  Eq. (1)  by  setting 


Y  X 

—       x  =  — 

L  '  L 


and  dividing  through  by  EI/L3  (where  L  =  length  of  the  beam) , 
we  obtain, 

2 
where  P  =  TL  /EI  =  nondimensional  force 

o      4 
Q  =  ft-Aoo  L  /EI  =  nondimensional  frequency. 

By  standard  methods,  we  obtain  a  general  solution 

y  =  C,Sinhax  +  C2Coshax  +  C  Singx  +  C.Cosftx  (2) 

i 

1  P2 
2   /  P 

(3) 


■i 


p2         p 

—  +  Q   -  — 
4*2 
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Constants,  C  ,  C  ,  C  ,  and  C.,  in  Equation  (2)  can  be 
evaluated  for  fixed-fixed,  free-free,  fixed-simple,  simple-simple, 
and  other  beam  end  support  conditions.   For  the  most  general  case, 
i.e.,  elastically  constrained  ends  condition,  no  results  are  found 
in  the  references.   This  case,  however,  is  the  most  important  for 
bridge  members,  and  we  shall  evaluate  the  constants  for  this  case 
presently. 

The  boundary  conditions  are 
y(0)  =  y(l)  =  0 

2 
d  y   _    dy 


i  -5lS  =0     at  x  =   0 


(4) 


dx        dx 


£3.   +  G  dy  =  o     at  x  =  1 
dx1  zdx 


kjL         k2L 

where     1~   EI   ,     2  = and  k,  ,  k„   are 

EI  12 


respective  torsional  moduli  at  the  ends  of  the  beam.   By 
applying  the  four  conditions  in  Eq. (4)  to  Eq.  (2) ,  the  vanish- 
ing coefficient  determinant  yields 

2     2  2  2    2 

{(a  +  6  )   +  0\02    (  a  -  3  )}   Sinha  Sm3  -  2  aia2a3 

2    2   r 
(Cosha  CosB  -  1)  +  (a i   +a2)  (a  +  3  )  {aCosha  Sin  3- 

BSinha  Cos3)  =  0 (5) 

Equation  (5)  is  the  most  general  frequency  equation  for 
small  amplitude  transverse  vibration  of  an  Euler-Bernoulli 
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beam.   As  far  as  we  know,  it  has  never  appeared  in  literature. 
We  can  show  that  it  includes  all  existing  results  as  special 
cases. 

After  the  natural  frequencies  are  evaluated  from  Eq. (5) , 
we  can  substitute  the  results  back  into  Eq. (4)  to  obtain  the 
relationship  between  the  constants  of  integration.   By  arbi- 
trarily setting  C„  =  1  and  knowing  C  =  -C„  we  can  determine 
uniquely  C.  and  CU.   Finally,  we  obtain  the  mode  shape  as 
follows:  (subscript  n  omitted) 
w(x)=  K  Sinh  ax  +  K2  (Cosh  ax  -   Cos  £x)  +  K3Sin  3x  ...(6) 


where 


K  = 

1 


K2  = 


K3  = 


2    2 
a  +  B 


ai3 


2         2  2 

-(a  Cosha  +3  Cos3)     ~3  Sin3+  C23Cos 


-  a2 (aSinha  +3Sin3) 


a^a 


ai3 


2  2 

a  Sinha  +  cuaCosha     -3  Sin3  +  a23Cos 


Oia. 


2    2 

a  +  3 


2  2         2 

a  Sinha  +  a2aCosha     -(a  Cosha  +3  Cos3) 


-a2(aSinha+  3Sin3) 

Equation  (6)  is  the  most  general  modal  shape  expression 
in  the  framwork  of  Euler-Bernoulli  theory.  Several  degenerated 
cases  are  interesting.   For  example: 
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For  a  beam  with  simple  supports  at  each  end, 

al  ~a2  =0 
so  Eq.  (5)  becomes 

2   2  ? 
(  a  +8  r-  Smha  SinB  =  0 

which  means   Sin  3=  0  ;   since  Sinha^  0 
therefore        8  =  n  tt  . 

The  associated  mode  shape  is  w  =  SinBx. 
For  the  case  of  a  beam  with  a  simple  support  at  one  end 
and  a  fixed  support  at  the  other  end, 


o1  -   0  ,      a2->  «> 
so  Eq.  (5)  becomes,  . 


ct 

_  tan 8=  tanha 


The  associated  mode  shape  is  w  =  8Cos8  Sinhax  -aCoshaSinBx 
For  a  beam  with  both  ends  fixed, 


so  Eq.  (5)  becomes, 


a    R 
SinhaSin8  -  — — — —  (CoshaCos8~  1)  =  0. 


The  associated  mode  shape  is 


w(x)  =  (aSinha  +8Sin8)-(8Sinhax  -  aSin8x)  -a8(Cosha-  Cos8)- 
(Coshax  -  Cos8x). 

Finally,  for  a  beam  simply-supported  on  one  end  and  partial- 
ly fixed  on  the  other ,  the  frequency  equation  becomes , 
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2    2 
(  a  +  3  )  Sinha  Sing  +a2(aCosha  Sing  -  gSinha  Cos  3}=  0 

and  the  corresponding  mode  shape  is 


2  2 

w  =  (  a28Cosg  -3  Sing)  Sinhax  -  (a  Sinha  +02aCosha) Singx 


we  have   computed   numerical  results  of  Eq. (5)  to  cover  all 
possible  conditions  of  beam  end  support,  for  cases  where  a,/  a9 
and  for  the  cases  of  applied  tension.   The  results  are  shown 
in  Appendix  B. 

(2)   Torsional  Vibrations 

Non-uniform  torsional  vibrations  of  a  thin,  open  sec- 
tion are  governed  by 

d4<S>       d2$ 

EC  — x  -   GC — 5"  +  Pi  W2$  =  0 
w  dX4       dX2      P 


where  EC  =  warping  rigidity 
w 

GC  =  torsional  rigidity 


I   =  polar  moment  of  inertia 
P 

$   =  angle  of  twist  per  unit  length, 


To  take  care  of  axial  tension  a  term 


T     d2$ 
A   p  dX2 


must  be  added  to  the  equation. 

$  X 

To  non-dimensionalize  it,  we  set  <f>= —  x  = — 

L  L 


dx4    1  dx2     1 
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where  P     =    (GC  +  TI  /A) L2/EC 
1  p  w 


Q,  pi    oj2L4/EC 

1  p  w 


The  general  solution  is  again  of  the  form  of  Eq.  (4) 
except  that ,  in  this  case , 


jp2  p 

'4     1      2 


(7) 


By  this  redefinition,  the  results  for  transverse  vi- 
bration can  be  re-interpreted  to  apply  also  to  the  correspond- 
ing cases  for  torsional  vibrations.   It  should  be  noted  that, 
a  simply  supported  end  in  torsion  means  no  restriction  to  warp- 
ing of  cross-section  at  that  end  while  a  fixed  end,  full-re- 
striction to  warping.   In  either  case,  the  angle  of  twist  is  zero. 

Some  contemporary  bridges  employ  HP  sections  as  slender 
members,  which  can  only  undergo  nonuniform  torsion  as  long  as 
the  joints  remain  unyielding  against  twisting.   Others  employ 
Box  Sections  which  can  resist  no  warping  and  wherein  the  uni- 
form torsion  prevails. 

It  is  interesting  to  compare  the  effect  of  warping  on 
the  natural  frequency  of  a  HP  beam  in  torsional  vibrations . 
For  sake  of  comparison,  we  consider  simply  supported  beams. 
Solving  f or w  in  Eq.  (7) 
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2    ■?  TI 

4    \r  n  tt^(GC  +  — rP-) 

2      n   tt  ECw  A 

w  =ptt4 — +  — zn? — 

P  p 


in  which  g  is  already  substituted  by  mr  .   Calculation  shows 
that  for  T  =  0,  the  fundamental  natural  frequencies  of  a  typical 
HP  section  can  possibly  be  raised  as  much  as  20%  due  to  warping. 

For  thin  closed  sections,  EC^  =  0,  and  the  differential 
equation  is  reduced  to  a  second-order  one.   Strictly  speaking, 
EC  is  not  zero,  but  is  very  small,  so  that  the  equation  is  a 
singular  fourth  order  one,  which  requires  a  singular  perturbation 
solution.   For  our  purpose,  we  can  ignore  the  "boundary  layer 
effect"   and  treat  it  as  a  second  order  equation.   It  can  be 
solved  in  a  straight  forward  manner,  but  we  omit  it  here,  since 
it  is  much  simpler  to  degenerate  the  frequency  equation.   See 
Part  C  ,  Chapter  II  for  details. 
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II.   RESPONSES  IN  SUBCRITICAL  FLOW 
(1)   Hartlen-Currie  Model 

The  response  of  an  elastic  beam  to  wind-induced  exci- 
tation should  be  determined  by  the  fluid-structure  interaction 
analysis  which  involves  the  integration  of  Navier-Stokes  equa- 
tion and  matching  of  pressure  and  displacement  interfaces .   A 
complete  analysis  is  beyond  the  state-of-the-art  or  rather, 
if  it  could  be  done  at  all,  the  expenditure  for  computer  usage 
would  be  prohibitive  for  the  present  purpose.   Therefore,  it  is 
expedient  to  resort  to  a  phenomenological  approach  as  described 
below. 

In  the  Strouhal  range,  the  amplitude  of  oscillation  of 
an  elastic  beam  is  nil  in  a  non-resonant  condition;  even  in  a 
resonant  condition,  i.e.,  the  Strouhal  frequency  is  close  to 
one  of  the  natural  frequencies  of  the  beam,  the  amplitude  re- 
mains quite  small,  on  the  order  of  ten  percent  of  the  dimension  'd'  of 
the  beam.   For  simplicity,  we  shall  neglect  the  variation  of 
the  response  of  the  beam  along  its  longitudinal  axis; a  segment  of  the 
beam  is  allowed  to  move  as  a  rigid  body  in  response  to  the  oscillat- 
ing lift  force.   The  elasticity  effect  is  represented  by  springs 
at  its  ends . 

Now  consider  an  elastically-supported,  viscously-damped, 
rigid  dxd  H-beam  of  length  L  (Fig.  31)  exposed  to  a  flow  of  free- 
stream  velocity  V.   The  beam  of  total  mass  m  is  elastically- 
supported  on  springs  of  stiffness  K  and  the  viscous*  damping  co- 


*  Hysteretic  or  structural  damping  introduces  only  minor  changes 
of  small  magnitude. 
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DETAIL  A  -  Hartlen-Currie  Model 
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Fig.  31    Vibrating  Beam  Modelled  by  Hartlen- 
Currie  Lumped  Mass  Theory 
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-efficient  C.   The  equation  of  motion  of  the  beam  is  given  by 

mw(t)   +   C  w(t)   +  K  w(t)   =  F(t)  (1) 

where  the  dot  over  the  displacement, w,  indicates  a  derivative 
with  respect  to  time  and  F  denotes  the  external  lift  force 
acting  on  the  beam  given  by 

F(t)  =  h   f^V2dLC    {2) 

wherein  C  denotes  a  time-dependent  lift  coefficient  and  pa  is  the  mass  density 
of  air. 

For  ease  of  subsequent  manipulations,  we  non-dimen- 

sionalize  Eqs .  (1)  and  (2)  by  substituting  Eq.  (2)  into  Eq.  (1) 

and  dividing  by  md  to  obtain 

2 
W"(t)   +  2YW'  (t)   +  w(t)  -  afi0-  C  (t) (3) 

J_i 

where  the  prime  indicates  differentiation  with  respect  to  the 
argument,  W  =  w/d,  %   =   /^}    T   =   ^  y=  Q/2m^    a  =  pd2L/87r2St2m> 

0,0  =   fs/fn=  St(V/£nd).   Such  that  u^  denotes  the  natural 
frequency  of  the  member  in  radians  per  second,  fs  denotes  the 

frequency  of  vortex  shedding  (f s  =  StV/d) ,  fn  is  the  natural 

th 

frequency  of  the  member  in  cycles  per  second  for  the  n   mode 

of  vibration,  and  St  denotes  the  Strouhal  number. 

The  nature  of  the  fluid  behavior  is  complex.   It  is 
generally  agreed,  however,  that  it  is  nonlinear  in  nature  as 
exemplified  by  the  "lock-in"  phenomenon.   We  recall  that  the  lock- 
in  phenomenon  describes  a  finite-spectrum  resonance  of  a  speci- 
men, in  contrast  to  a  point- spec trum  or  a  delta  function  (spikes) 
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type  of  resonance  which  is  characteristic  for  a  linear  system. 

In  a  linear  system,  the  resonance  will  occur  only  if 
the  Strouhal  frequency  is  exactly  equal  to  the  natural  frequency 
of  the  elastic  system.   In  a  nonlinear  system,  on  the  other  hand, 
the  resonance  will  begin  even  when  the  Strouhal  frequency  is 
somewhat  less  than  the  natural  frequency  and  will  persist  even 
after  the  Strouhal  frequency  exceeds  the  natural  frequency  by 
a  substantial  margin.   Within  this  range,  the  natural  and  the 
Strouhal  frequencies  lock-in  to  a  common  frequency ,  which  is 
close  to  the  natural  frequency  of  the  system. 

Following  Hartlen-Currie,  we  postulate  that  the  fluid 
action  obeys  the  best-known  and  perhaps  the  simplest  nonlinear 
oscillator,  the  van  der  Pol  model: 

c"  (t)  -efi0C'  (r)  +  -  [  c;  (t)]  3  +fincT  (t)  =bW'(T)  ...(4) 
L       L      ao 

The  van  der.  Pol  parameters  e  and  £  ,  and  the  coupling  coefficient 
b  are  assumed  to  be  constants. 

(2)   Solution  of  Hartlen-Currie  Equations 

By  assuming  that  the  force  excitation  and  the  beam  dis- 
placement response  oscillate  at  the  same  circular  frequency   w 

c 

(close  to  co  )  ,  that  both  have  near  sinusoidal  waveform  with  little 
n 

amplitude  modulation,  and  that  the  excitation  leads  the  response 
by  some  phase  angle  <f>,  we  seek  solutions  of  Eqs.  (3)  and  (4)  of 
the  form 
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W  =  W  Sin  fix         (5a) 

C_  =  rsin(fix  +  $)    (5b) 

Lj 

fi  =0)  /W  fRrO 

c  n  VJl-' 


where,  to  the  zeroth  order  of  approximation,  w  an<^  T  are  constant 
amplitudes . 

Using  Eqs.  (5  a»b,e)  in  Eq.  (3)  and  equating  coefficients  appro- 
priately shows  that  C  and  T   can  be  determined  in  terms  of  the  as 

L 

yet   undertermined   Q,  and  W     through 

2    _           2 
TCos    fix   =    (1-    0,   )W/   afi0       (6a) 

2 
TSin  fix  =  2  /UQW/  aftQ       (6b) 

Substituting  Eqs (6a,£)into  Eq.  (4)  and  equating  the  like  co- 
effecients  of  the  trigonometric  functions  we  obtain, 

~,  2     2V  2       2     3CW2   3       2 

2Y«.-(n0  -  n  )-  en0n(i  -  n  )  -abfi0«  +-~9-5  n  (i  -  n  )z  =  o  ...(7a) 

4a  fiQ 
p    2    2  2   3CW2      t+ 

a  -  n)  (n0  -n  )  +  2  eYfiofi  " — 2~s  2y^  z  =0   (7b) 

where  Z  =  (1  -  ft2)  2+  4  y2  ft2 . 

~2 
Eliminating  W   and  understanding  that  Z^O,  we  obtain 

2            2                            2    2 
(fiQ   -   fi   )Z-   2abyfi0fi-  -0      ....(8a) 

2     _L 

or       Q0/fi  =  (1  -  2abYfi  /Z)  2 (8b) 

We  note  in  passing  that  there  are  three  constants  to  be 
calibrated:  b,  £  and  C.   The  additional  condition  can  be  obtained 
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from  the  complementary  solution  of  Eq.  (4) .   By  setting  the 
left-hand  side  of  Eq. (4)  equal  to  zero,  it  is  well-known  that 
the  solution  is 


-    3  e 

C  =1 (9) 

L  \/4  ? 


which  provides  the  needed  additional  condition  to  evaluate  the 
three  unknown  constants. 

(3)   Calibration  to  Narita  Data 

The  original  Hartlen-Currie  model  was  calibrated  for  a 
vibrating  cylinder  with  satisfactory  results.   We  propose  to 
calibrate  the  model  for  a  vibrating  H-beam.   The  only  known  H  beam 
experiment  applicable  to  the  present  case  was  conducted  by  Narita. 
We  shall  summarize  the  finding  as  follows : 

p32L  _o 

y   = =   5.59  x  10  J 

M       m 

y    =   2.67  x  10~2 


CT   =   1.15 


w   /d*  =   0.09   @  Qq=    1.06 
max 


From  the  above  data  we  can  compute 


V  7.08     ..-5 

a  =  , — - - =        x  10 


8  IT2  St2     St2" 


*Note:   The  deflection  arrived  at  using  the  Hartlen  Currie  Model  represents  the 

mean  deflection  of  the  beam.   So  W_av  should  be  considered  as  W    as  shown 
.   _ .    , ,  max  max 

in  Fig.  11. 
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Narita  did  not  calculate  the  Strouhal  number  in  this  case  which  we 

estimated  to  be  in  the  neighborhood  of  0.125  to  0.134.   Through  a  parametric 

study  the  best  fit  for  Eqs. (7b)  and  (9)  was  given  by 

St  =  0.130 
b  =  1.5 
£  =  0.5 

Figures  32,33,34  show  how  these  data  correspond  to  Narita' s  experimental  results. 

Figure  33  shows  that  with  b=7.5  there  was  excellent  agreement  between 

experimental  observations  and  analytical  predictions  of  the  frequency  data. 

For  the  more  important  amplitude  response  analysis,  however,  it  was  found 

that  a  value  of  b=l„5  and  £=0.5  provided  the  closest  fit  to  Narita' s 

observed  maximum  amplitude  response  of  0.0925  at  Q,Q=1.06.      Figure  34  shows 

the  result  of  the  parametric  study  for  frequency  response  using  b=1.5 

while  Figure  31  show  the  amplitude  response  correlation  between  experimental 

and  analytical  results.   The  analytical  results  note  a  maximum  amplitude 

of  0.1092  at  fiQ  =  1.06  which  will  yield  slightly  higher  deflections  (on 

the  order  of  18%)  over  Narita' s  experimental  observations.   This  will 

provide  an  added  safety  factor  in  analyzing  and  designing  for  maximum 

deflection. 

(4)   Comments 

It  is  seen  that  Hartlen-Currie  Model  also  gives  a  satisfactory  result 

for  vibrating  H  beam.   Its  accuracy  leaves  much  room  for  improvement.   One 

possible  candidate  is  Skop-Griffin  Model  which  contains  two  extra  parameters. 

This  model  requires  much  more  labor  in  calibration.   Other  shortcomings  of 

Hartlen-Currie  Model  are: 

(1)  it  predicts  a  frequency  hysteresis  loop  which  has  never  been 
observed  in  the  experiments 

(2)  it  fails  to  yield  accurate  information  beyond  the  maximum 
amplitude;  in  fact,  W  will  eventually  become  an  imaginary 
number  which  is  absurd. 

Figure  11  shows  a  relation  between  the  amplitude  and  the  damping 

coefficient  which  may  be  useful  in  studying  the  real  beam  response. 
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III.   RESPONSES  IN  SUPERCRITICAL  FLOW 

The  equation  of  motion  of  a  continuous  system  takes  the 
general  form 

m(X)w(X,t)  +  c(X)w(X,t)  +  L(w(X,t))  =  F(X,t)    (1) 

where  X  is  the  coordinate  position  of  a  generic  point  in  the 
continuous  structure,  m(X)  is  the  mass  density  per  unit  length, 
w(X,t)  is  the  displacement,  c(X)  is  a  viscous  damping  coefficient  , 
L(°)  is  a  self-adjoint  linear  differential  operator,  and  F(X,t) 
is  an  arbitrary  forcing  function.   Note  that  X  and  w(X,t) 
are  one  dimensional  in  the  present  context  but  may  be  inter- 
preted as  vectors  in  two  or  three  dimension  cases. 

Examples  of  the  form  that  L  may  take  are 

2 
3 
L  ( "  )  =  T— r(')  for  a  vibrating  string, 

dx2 

2       2 

3      3 


(EI— r('))        for  a  vibrating  beam, 


3X^    3X 

2      2  ? 

3     3  3  " 

and  =  ——EI — r(  * )  -  T  — ~~  for  a  vibrating  beam 

3X2   3X2        3X2 
with  tension. 

To  solve  (1)  we  assume  sufficient  conditions  to  insure 

w(X,t)  has  an  orthocronal  expansion  in  normal  modes  w  (X)  of 

"  n 

oo 

the  form  w(X,t)  =  I'  w  (X)q  (t)  ,  where  the  generalized  coordinates 

n=l 
q  (t)  remain  to  be   determined.   Using  this  expansion  we  may 


JL 

Hysteretic  or  structural  damping  introduces  only  minor  changes. 


00 


rewrite  (1)  as 


m(X)Z  w  (X)q  (t)  +  c(X)E  w  (X)q(t)  +  E  L (w(X) )q„ (t) 
n=l  n    n         n=l  n         n=l   n     n 

=F(X,t).  (2) 


Pick  a  particular  normal  mode  w  (X) ,  multiply  (2)  by  it 
and  integrate  with  respect  to  X  to  obtain  the  infinite  system 
of  ordinary  differential  equations  in  t, 


2.2 

_ .     f   < 
r   r  r   r   r  Jr  r  r  ar 


M   q      +    2ttM    f    y    q      +   M    4t    f    q      =    P     (t)  ,  (r=l,2...)     (3) 

•   r  r    r   r   r  v-  r^r  r 


where  M      =    X£m(X)w2(X)    dX, 

r  0  j; 

2ttM    f    y^   =    /£c(X)w2(X)     dX, 
r   r    r  o  r 

and  P    (t)    =    /   F(X,t)w    (X)    dX. 

r  °  r 

We  point  out  that  to  obtain  (3) ,  we  assumed  the  orthogonality 
conditions 

/J"m(X)w  (X)w  (X)  dX  =  5   M 
u     r    n  rn  r 

£ 
and  i\c(X)w    (X)w    (X)    dX  =    <5    4y    ttM   f         (i.e.  .we    assume 

0  r  n  rri     r       r   r  ' 

no  frictional  coupling  between  the  modes)   where  y   is  the 

J. 

generalized  damping  term,  commonly  known  as  the "viscous  damping" 
(i.e.,  it  radians  out-of -phase  with  the  velocity  w  but,  unlike 
the  hydrodynamic  damping,  is  independent  of  the  magnitude  of 
the  velocity  and  the  f  ' s  are  the  natural  frequencies  derived 

fromtherelationshipL(x)    =  .,2rf    )  \   . 

^  n         n    n 

Equations  (3)  govern  each  individual  natural  mode  of 
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vibration  and  consequently,  we  have  in  effect  uncoupled  all 
the  modes  of  vibration.   By  design,  equations  (3)  have  the 
same  form  as  equations  of  motion  of  a  simple  oscillator.   Thus 
methods  which  have  been  previously  developed  for  the  treatment 
of  random  vibrations  of  discrete  mechanical  systems  can  be 
used  here  in  developing  analogous  techniques  for  the  treatment 
of  random  vibrations  of  continuous  systems. 

Returning  to  the  solving  of  system  (3) ,  take  the  Fourier 
transform  and  obtain 


g    (f)    =   H(f)Pr(f) 

-i2frft 

—  oo 

where  qr(f)    =   L^^^  dt  > 

-i2irft 
Vf)    =   Z^Pr(t)e  dt, 


and  H(f) 


4TT2Mr.{(f2    -    f2)+   iy    f    f } 
r        r 

i  2  it  ft 


r    r  L    r 


Let  h(t)  =  /0H(f)e      df.   Then  we  obtain  the  solution 

q  (t)  =  /th(t  -  t)P  (t)  dx 
Jr      —  °°         r 

00 

=  /  h(x,)P  (t  -  t,)  dr. •  (4) 

0     1   r        l  1 


The  next  step  in  the  treatment  of  this  continuous  system 
is  the  development  of  relationships  among  the  statistical 
properties  of  the  physical  and  modal   excitations.   These 
are  of  primary  concern  because  the  original  influences  on  the 
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structure  were  not  specified  in  any  manner  and  so  are  assumed 
to  be  random  processes.   To  simplify  the  considerations,  we 
assume  these  random  processes  to  be  stationary,  i.e.  invariant 
under  time  shifts.   It  follows  then  that  the  autocorrelation 
function  F  (X  ,t)  associated  with  the  solution  of  (1)  will 
be  a  function  of  only  the  time  lag  x  for  any  given  point  X,. 
Specifically,  we  have  (where  E(')  represents  expected  values) 
that 

Rw(Xa,t)  =  E{w(XA,t)   w(XA,t  +  t)} 

=  Z  Ew„(X  )w  (X  )   R     (t)  (5) 

r  s  r   A   s   A    q,rs 

where       R     (t)  =  E{q  (t)q  (t  +  x)}- 
q,rs  r     s 

We  must  calculate  R     (x) .   To  wit,  since 

q,rs 

q    (t)q    (t   +    x)    =    l"    /"h'(T.)h(T1)P    (t    -    TjP_(t    - 

r  s  oolir  is 

x2    +    x)    dx    dx2 

where  P    (t   -    x)P    (t   -    x0    +    x)    =    /      f„F(S    'fc   _tJw    <xi  ) 

r  l      s  2  °°x  ir1 

F(X„,t    -    t,    +    x)w    (X    )      dXn   dX 
2  1  s       2  1         2 

then  R     (x)  =  /"  /"h(T1)h(T0){/J  fZR(XlfX?; 

q ,rs  oo  L  2         oofi^ 


x)w    (XjwJXJdX  dXJdT.dTj 
rls^        1  lz 


where  R  (X  ,X  ;x)  is  the  autocorrelation  function  associated 
p   1   2 

with  the  forcing  term  F(X,t). 

Next  we  wish  to  determine  the  spectral  density  function 
S  (X  ;f)  which  corresponds  to  the  autocorrelation  function  Rw(X1,x) 
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This  is  done  by  taking  the  Fourier  transform  of  (5)  and  using 

i2frf  x 

00 

the  Wiener-Khinchin  relation   S(f)  =  /  2e     R(x)dx   between 

—  00 

a  spectral  density  S(f)  and  its  associated  autocorrelation 
function  R(t)  . 
We  get 

S  (X  ;f)  =  I    Zw  (Xn)w  (X  )S     (f)  (6) 

wl       rsrlsi   q,rs 

where  S^   ^(f)    =   H    (-if)    H    (if)     {/°°    /~S„ (X, ,X_ ; t) w 

q,rs  r  s  ooF12  r 

(xx)w   (x2)   dx1dx2.} 

S  (X  ,X9;t)  is  naturally,  the  spectral  density  associated  with 

F   1   ^ 

Rp  (X.. ,X„ ; x )  . 

In  our  case,  the  frequencies  are  distinct  and  the  damping 
is  small  so  that  the  peaks  are  well  separated.   We  can  approximate 
by  neglecting  the  cross  products  and  obtain  the  simpler  expression 

S(X,,f)  =  Ew2(X  )S    (f) 

w   A      r  r   A  q  r 

where  S    (f)  =  |H  (f )  |   *  {/*  /Jwr  (Xj)wr  (X2)  Sp  (X^X^-f )  dX^Xj  }  • 

If  we  assume  that  the  lift  force  is  spatially  independent 
i.e.  if  there  is  perfect  coherence  spanwise,  then  S  (X,X2;f)  = 
SF(f)  and  hence 

S    (f)  =   H  (f)       {/£  .Tw  (X-,)w   (X9)S  (f)dX-.dX2}  • 
q,r       '  r  o   o  r   l   r    2   p      1^ 

It  therefore  follows (by  taking  a  mean  value)  that 

E{q2(t)}  =  /°°S  (f)    (/Jw  (X)dX)2  *|H  (f)  I   df 
r         o  F        u  r  r 

=  Sp(f)|     (/£w  (X)dX)2/ro|H  (f)|2df. 

r     f=f    °  r         or 
^  -"-o 
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00  r\ 

We  must  now  evaluate    /  |h  (f) I  df.   Using  the  method  of 

o  '  r    ' 

residues,  one  finally  gets 

2  1 

/£|H  (f) |   df  = 


64^r  MrYrfr 


(In  the  engineering  analysis  the  hysteretic  damping 

or  structural  damping  is  denoted  by   n   which  should  satisfy 

r 

2 

M  qv  +  Mr(2-rrfr)  {1  +  in  }qr.    =   P„(t).  ._  . 

r  r    r    r         r  ^r     r  (3a) 

In  this  case 


H(f)  = 


4tt  M  {  (f2  -  f2)  +  in  f 2 
r    r  r  r 

2  1 

so  that     |H(f)  |   =     2    o    o  9  ?     24 

(4tt  M  )z{  (fz   -  fz)      +   n  f  > 
r     r  £   r 


which  for   nr.<<l  reduces  approximately  to 


r 

1 


/°°|H(f)  |2df  = 


o1  M  'i  7  1 

4M^(2irf   Jn 
r    r    r 

which  correlates  nicely  with  the  preceding.)   Now  we  invoke 

the  fact  that  for  sharp-peaked  |H(f) I2,  f   is  equal  to  f   approx- 

o  r   ^ 

imately  and  that  the  rrns  qr(t)  is  the  square  root  of  E{q£(t)}. 

Thus  we  have  as  an  expression  for  the  root  mean  square  of  q    (t)  '. 


/*w  (x)  dx  /s~7fTJ 
rms  q  (t)  — 


r      J  T[°   I"   Y 

r   r 


-105- 


Since  S„(f  )  is  usually  normalized,  i.e. 

r    r 

/"Sf (fr)  df  =  /~Sf  (st)d(st)  =  1   where  St  is  the 
Strouhal  number  and 

SF(fr)  =  tepCLdV2)2  Sf(fr)  ,     we  obtain 


•L, 


rms(q    )    =      2  ihp£i$VZ) 


8Mr  ^  7rfr)  3Yr 


Finally,  returning  to  the  solution  S(X,t), 


rms(w)|  =    Zw    (X    )     *    rms (q    ) 

A 


L 


r      A         or  ,,__„    JTTzv     /    f      r      r 


=    r"T"  3/2"    2        (^CLdV    >w- 

r    8M     (     tt)-3/^     fj  V  y 


■/ 


r 


wr(XA)     /0w       dX                                   /S     (St)St 
j  -7       3/2       2     (PaCLdv2)    V~~ (?) 

i6Mrc  tt)      frz  y        Yr 


rms   M 


2 
8  w 


=    rms    EI 


A  d  x 


X — Xa 


(8) 


To  evaluate  Equation  (8) ,  the  mode  shapes  must  be  twice 
differentiated,  which  is  complicated  to  evaluate  numerically, 
except  for  simply-supported  beams.   In  general,  it  is  simpler  to 
use  the  deflection-moment  relationship  as  derived  in  a  later 
section.   (See  Part  C  ,  Chapter  III.) 


r  9 

The  quantity  (^pCTdV  )  is  the  lifting  force. 

a  A4 
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PART  C 
DESIGN  PROCEDURES  AND  NUMERICAL  EXAMPLES 

I.  NATURAL  FREQUENCIES,  MODES  AND  WIND  VELOCITIES 

{ 1 )  General  Case  for  Bending  and  Torsion  for  Thin  Open  Sections 
The  general  equations  developed  above  were  used  to  develop  the 
following  design  procedure.  The  designer,  having  made  a  preliminary  choice 
of  a  beam  for  use  in  the  structure,  can  easily  determine  the  wind  velo- 
cities which  will  cause  the  beam  to  vibrate  in  resonance.   First,  the 
value  of  v/p  for  bending  can  be  calculated  as : 

for  open  or  closed  cross  sections 
where     T  =  Tensile  force  on  beam 
L  =  Length  of  beam 
E  =  Modulus  of  elasticity 
I  =  Moment  of  inertia  about  axis  of  vibration 

and  for  torsion: , 

f[QC   +(T/A)lp]L2 

/P  =/ — ' — — — -   for  open  cross  sections. 

<  ECw 

Showing  the  value  of  P  the  designer  can  go  to  the  provided  graphs  of 

VP~ vs.  vQ  (Figs.  1-10)   and  determine  the  values  of  /Q  for  the  first  three 

modes  of  vibration  on  the  given  beam  of  open  cross  section.   These  graphs* 

consider  the  end  support  conditions  of  the  beam  and  provide  results  for 

beams  with  both  ends  having  simple  supports,  both  ends  having  fixed 

supports,  and  a  wide  variety  of  end  support  combinations  between  these 

two  extreme  conditions. For  beams  of  closed  cross  section,  the  designer 

can  use  the  average  values  of  /Q  given  in  Table  1. 


*  For  more  complete  graphs  including  unsymmetric  cases,  see  Appendix  B, 
Figs.  36  -  84. 


-107- 


The  value  of  /Q  is  related  to  the  natural  frequency  of  the  beam 
as  follows: 


PmA-nL4 


for  bendinq  of  all  cross  sections 


EI         (open  or  closed) 


where 


pm=  Mass  density 

A  =  Cross  sectional  area 

(jo  =  Natural  frequency 
n 

L  =  Length  of  beam 

E  =  Modulus  of  elasticity 

I  =  Moment  of  inertia  of  section  about  axis  of  vibration 

n  =  Mode  of  vibration     (n  =  1,2,3) 


For  torsion  of  open  cross  section 


or  /Q  = / for  torsion  of  closed  cross  section- 

GC 
(For  derivation  see  next  section.) 

Where : 

I   =  Polar  moment  of  inertia 
P 

C  =  Warping  constant 
G  =  Modulus  of  rigidity 
C  =  Torsional  constant 
For  the  beam  with  both  ends  simply  supported  the  relationship 

between  ,/P  and  /Q  is  simple  and  is  given  as; 

Ah  n  2 

Q  =  n*  ir  +  nziT  P         n  =  1,2,3 
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However  the  relationship  between  P  &  Q  becomes  complicated  for 
other  end  support  conditions.   Therefore,  the  graphs  (Fig.  1  -  10)  have 
been  worked  out  by  computer  analysis  to  provide  the  designer  with  a 
simple  means  of  finding  the  value  of  /Q. 

After  determining  the  values  of  /Q  for  the  first  three  modes 

of  vibration  from  the  graphs  the  designer  can  use  them  to  directly 

determine  the  natural  frequencies  that  correspond  to  them.   The 

designer  can  find  the  Strouhal  number  corresponding  to  the  shape  of 

his  beam  from  the  material  provided  in  Chapter  IV,  part  A .   The  wind 

velocity  related  to  the  natural  frequency  and  Strouhal  number  is: 

0)  d 
V=  -£- 


2  TT  St 


Where: 


V  =  Wind  velocity 
St  =  Strouhal  number 


a)  =  Natural  frequency 


d  =   Depth  of  section  transverse  to  wind  direction 
Alternately,  the  designer  can  compute  the  wind  velocity  directly, 
knowing  the  Strouhal  number  and  the  values  of  /Q.   He  must  determine 
a  value  of  Z  which  relates  /Q  and  St  as  follows: 

For  bending  of  open  or  closed  cross  sections,  we  know  that: 

2  4 

pAto  L  w  d 

Q  =  ,   and    St  =  

EI  2ttV 


/^A1    L  2ttV 

therefore ,       /q  -  I  —   — St 

/  EI      d 
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or,  we  define      Z  =- 


v/st 


4  , 2~> 

'PAL4V24tt 
Z  =" 


Eld2 


A  similar  approach  can  be  taken  to  show  that  for  torsion: 
4 


1  2041 

pi    4it    \TL 
mp 


For  open  cross  sections 


or 


/pi    4tt   V2L2 
mP 


GCd2 


For  thin  closed  cross  sections  , 


Computing  the  values  Z  =  /§/  /St/   or  from  Fig.  35  ,  and  using 
the  definition  of  Z,  the  designer  can  easily  find  the  wind  velocities 
corresponding  to  each  mode  of  vibration  from  the  following  relations: 

2    d    P1 
V  =  Z   /  —     For  bending  of  open  or  closed  cross 

2ttL  7  nA      sections 
m 

and 


-\ 


J  EC 
2         /   w 
V  =  Z   — —  / For  torsion  of  open  cross  sections 


2 

2ttL   /  pi 


m 


P 


or 


For  torsion  of  closed  cross  sections. 


Knowing  the  wind  velocity  which  will  cause  vibration  in  the  beam 
being  considered  and  the  weather  data  of  wind  velocities  in  the  de- 
sign area,   the  designer  can  decide  whether  the  selected  beam  is 

sufficient  or  whether  another  section  should  be  used. 

Examples  of  this  design  procedure  can  be  found  in  Examples  1  to  4, 
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(2)  Natural  Frequencies  for  Torsional  Vibration  of  Thin  Closed  Section 

The  governing  equation  for  bending  of  a  structural  member 
with  a  thin,  closed  cross  section  such  as  a  box  or  circular  pipe 
section,  is  the  same  as  that  for  a  member  with  an  open  section  and 
the  value  of  /Q   can  be  obtained  as  usual  from  the  graphs  provided 
(Figs. 1-10) .   However,  for  a  member  with  a  thin,  closed  cross  section 
in  torsion/ C  equals  zero.   The  general  governing  equation: 

PlfoA4  ._  4  4      2  2   [GC  +(T/A)IplL2 


=n  tt  +  n  t 


EC  EC 

w  w 


reduces  to : 


pi  (/  L2         o      o   (T/A)  I 
nip           2   2     9  2        P 
c =  n^  tt  +  n"S   — 


GC  GC 

which  is  of  the  form  Q  =  n27r  +  n2-rr  P.   Therefore,  for  thin,  closed 
sections,  these  parameters  are  non-dimensionalized, 

(T/A)  I 


Q  = 


GC 


2  2 
pi  OL,  IT 
tn  p  n 


GC 

With  this  form  of  P,  the  values  of  /P     can  easily  be  calculated 
for  any  member  being  analyzed  just  as  in  the  case  of  an  open  cross 
section. 

Within  the  range  of  allowable  tensile  stress  (from  T/A  =  0  ksi 
to  T/A  =  60  ksi)*,  the  values  of  /Q  for  thin,  closed  cross  sections 


*  American  Institute  of  Steel  Construction  -  Manual  of  Steel  Con- 
struction. 
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vary  over  a  small  range.   Over  this  small  range  the  average  values 
of  /Q     are : 

Mode  1,  fQ   =  1.77 

Mode  2,   /Q  =  2.51 

Mode  3,   /Q  =   3.07 
which  correspond  to  a  tensile  stress  of  T/A  =  20  ksi. 

These  values  can  be  used  for  all  shapes  of  thin, closed  cross 
sections  in  torsion  having  any  specified  axial  tension  with  an 
expected  error  of  well  below  1.0%. 
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Fig,  35  Velocity  Coefficients  Z  Related  to  Q  8  St 
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(3)  DESIGN  EXAMPLES 

This  section  will  provide  practical  examples  that  demonstrate 
the  application  of  the  previously  presented  theory.   The  examples 
assume  that  the  designer  has  already  chosen  a  particular  structural 
shape  for  use  and  wishes  to  calculate  the  wind  velocities  which 
correspond  to  the  member's  first  three  modes  of  vibration. 

Depending  on  how  the  structural  members  are  connected,  the 
designer  has  to  determine  the  end  support  conditions  of  the  member.   After 
deciding  this,  simple  calculation  and  the  charts  provided  within  this 
report  will  quickly  lead  to  the  determination  of  the  wind  velocities  that 
will  directly  affect  the  member. 

Historical  data  of  wind  velocities  in  the  vicinity  of  the  structure 
can  be  obtained  and  consulted  giving  the  designer  an  accurate  under- 
standing of  whether  the  given  member  will  be  affected  by  existing  site 
winds.   If  the  wind  velocities  commonly  experienced  in  the  site  area 
correspond  with  those  calculated  the  designer  may  decide  either  to 
choose  a  different  member  for  the  structure  or  to  modify  the  structure 
using  the  methods  discussed  within  this  report  to  insure  minimum  wind 
induced  vibration. 

Previous  studies  have  shown  that  circular  cylinders  behave  in  an 
unusual  manner  once  the  Reynolds  number  of  the  airstream  flow  reaches  a 
value  above  1.0  x  10  .  At  this  point  the  member  enters  a  transitional 
range  which  extends  to  a  Reynolds  number  of  about  4  x  10^  depending  on 


.m  . 


the  upstream  turbulence  and  local  disturbances.   Within  this  range, 
the  circular  cylinder  will  jump  from  a  subcritical  state  where  vibrations 
concentrate  at  the  fundamental  natural  frequency  of  the  member  to  a 
supercritical  state  where  random  vibration  occurs.   In  the  supercritical 
state  vibrations  have  wide  dispersion  and  a  wind  velocity  corresponding 
to  the  fundamental  mode  of  vibration  excites  the  third  mode  of  vibration 
of  the  member. 

Whether  such  an  occurance  is  possible  for  structural  cross  sections 
other  than  circular  was  not  evident  in  the  literature  researched  for 
this  report.   If  the  future  research  should  show  this  phenomenon  to  be 
a  general  rule  then  the  procedure  outlined  in  example  5  can  be  applied 
for  checking  the  member  in  the  supercritical  range.   Otherwise,  the 
subcritcal  flow  will  be  the  dominant  criterion  of  analysis  regardless  of 
the  Reynolds  numbers.   After  calculating  the  wind  velocities  that  will 
directly  affect  the  chosen  member,  the  designer  can  determine  whether 
the  member  will  undergo  subcritical  or  supercritical  vibration.   If  the 
flow  is  in  the  transitional  or  supercritical  range  the  member  should 
be  analyzed  considering  random  vibrational  effects  with  modification 
made  accordingly. 

It  is  therefore,  the  aim  of  these  examples  to  provide  insight  on 
how  problems  can  be  realistically  evaluated.  Various  end  conditions, 
open  and  closed  cross  sections,  and  vibration  levels  will  be  examined. 

Hopefully  these  examples  will  provide  significant  demonstration  of 
the  theory  developed  earlier  and  its  application  in  the  design  process 
to  reduce  or  avoid  possible  dangerous  affects  of  wind  induced  vibrations. 
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DESIGN  EXAMPLE  NO.  1 

Assume  a  beam  HP  12  x  74  is  to  be  examined.   It  is  45  feet  long, 
carries  a  tensile  load  of  87.2  KIPS,  and  has  simple  supports  at  each 
end.   (0|  =  (72  =  0) 

With  this  preliminary  information,  the  designer  can  begin  the 

analysis.   First  the  properties  of  the  HP  12  x  74  section  should  be 

determined.   The  Manual  of  Steel  Construction  put  out  by  the  American 

Institute  of  Steel  Construction  provides  easy  access  to  such  information 
Y 


„tf 


Properties  for  Designing  (AISC  Manual  p. 1-48) 

HP  12  x  74  section 

2 
Area  =  21.8  in. 

Depth  (D)  =  12.12  in. 

Width  (bj  =  12.217  in. 
f 

Thickness : 

Flange  (t  )  =  0.607  in. 

Web  (t  )  =  0.607  in. 
w 


Orientation  in  the  wind 


Elastic  Properties: 


k 


h=D 


U_l 


xx  -  axis 


yy  -  axis 


d=b, 


I  =  566  in. 
S  =  93.4  in. 
r  =  5.10  in. 


I  =  185  in. ' 


s  = 

30.2   in. 

r  = 

2.91   in. 

I  = 
P 

751   in.4 

Modulus  of  elasticity  (E)  =  29  x  10  ksi 
Warping  constant (Cw)  =  6120  in. 6 
Torsional  constant  (C)  =  2.97  in.4 
Poisson's  ratio  (r)  =  0.3  (for  steel) 
Modulus  of  rigidity  (G)  =  E 


1.115  x  1(T  ksi 


2(l+r) 
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Once  the  specific  properties  for  the  section  are  available,  the 
designer  can  begin  the  calcu  lations  by  solving  for  the  bending  and 
torsional  values  of  /?. 

For  Bending:  2       2 

TL2  =   (87.2  kips) (45ft.)  (144in'/ft2)   =  4  74 
P  = 


Elyy*    (29  x  103ksi) (185  in.4) 


v^"D  =  2.177 

D 


for  torsion 

2 


(GC  +  IpT/A)L' 


PT  = 


Ff 
c^w 


87.2  kips        ?    .  2 

(1.115X  104ksi)(2.97  in.4)  +  (751  in.4)( )  (45ft)^(144in/  2) 

21.8  in.2  " 


P  =  5.936  x  10 


»/p  =  7.705 


(29  x  103ksi)(6120  in.6) 
1 


With  the  above  values  of  /P  for  bending  and  torsion  and  with  knowledge 
of  the  end  support  conditions,  the  designer  can  easily  read  the  values  of 

A  ** 

vQ  for  the  first  three  modes  of  vibration  from  Figure  1,  page  15. 


Vibrations  will  occur  about  yy  -  axis  for  orientation  of  member  in  this 
example. 

•k  -k 

Note:   For  cases  where  the  member  is  simply  supported  at  both  ends, 
the  relationship  between  /F  §  vQ  reduces  to: 

^Q  =  nV  +  n2Tr2  /P 

If  desired,  the  designer  can  directly  calculate  vQ  but  the  charts  provided 
are  accurate  and  easily  used. 
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MODE 
n 

BENDING 

TORSION 

1q 

1 

2 
3 

3.49 
6.49 
9.55 

5.00 

7.90 

10.70 

The  next  step  is  to  calculate  a  value  Z     for  both  bending  and  torsion 

This  value  will  be  used  to  calculate  the  velocity  of  wind  correspondin 

the  first  three  modes  of  vibration.   The  designer  must  choose  an 

appropriate  Strouhal  number  for  the  section  being  analyzed.   For  this 

example  the  Strouhal  number  for  a  "H"  shaped  cross  section  is: 

S  =  0.12* 
t 

** 
then    Z  can  be     calculated       from  the  following  relationship: 

Z    =   ^Q  ^    v/s^ 

and  the  values  of  Z  for  this  example  are  as  follows: 


MODE 
n 

BENDING 
ZB 

TORSION 
ZT 

1 

2 
3 

10.075 
18.735 
27.569 

14.434 
22.805 
30.888 

Now  the  wind  velocities  can  be  computed  using  the  relationships 
For  bending: 


Reference  material  within  this  report  provides  Strouhal  number  information 
for  various  cross  sections. 

** 

Values  of  Z  can  also  be  read  directly  from  chart  provided  (Fig.  35, 

p.  113)  if  desired. 
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where:   I  =  The  moment  of  inertia  along  the  axis  about  which 
the  vibration  takes  place. (For  this  example  and 
section  orientation,  use  Iyy) 

o  =  Mass  density  of  the  material  of  the  member.   For 
this  example  the  material  is  steel  with  a  density 
of  489  / n.T,,     the  mass  density  is  therefore 

„„  lb-sec 


ft' 


V 

z2 

T  " 

d 

4891b/ 
P      =               Xft3              =      1! 
m                    tt 

32  ft/        2 
JZ     /sec 

2 

p      =   7.324  x   10~      lb_sec 

4H                                                               ... 

in.4 

torsion: 

1    ECW     ' 

2ttI/ 


tfi 


where:   I  =  The  polar  moment  of  inertia  of  the  section  (I  =  I   +  I   ) 
P  p    xx    yy 

p  =  Mass  density  of  the  section 

In  J 

The  values  of  V   for  this  example  have  been  carefully  calculated  and 
are  summarized  in  the  following  table.   Along  with  the  Reynolds  numbers 
that  correspond  to  each  velocity.  The  value  of  the  Reynolds  numbers 
are  arrived  at  through  the  simple  calculation  of: 

Re  -^— 

en     v 


where  u  =  Kinematic  viscosity  of  air 

sec- 


(u  =  2.3  x  10"2  in. 2/   ) 
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MODE 

BENDING 

TORSION 

n 

Wind 
Velocity 

V 

mph (km/hr) 

Reynolds 
Number 
Re 

Wind 
Velocity 

v 

mph (km/hr) 

Reynolds 

Number 
Re 

1 

2 
3 

22.292(35.876) 
77.088(124.06) 
166.92(286.63) 

2.084xl05 
7.207x10s 
15.605x10s 

44.837(72.16) 

111.93(180.14) 

205.336(330.46) 

4.19xl05 
10.464x10s 
19.196xl0b 

With  the  above  information  relating  wind  velocity  to  the  natural 
frequency  of  the  member,  the  designer  is  armed  with  sufficient  information  to 
make  a  decision  as  to  the  choice  of  this  member  for  the  structural  examining 
of  the  resulting  wind  velocities.  The  designer  can  readily  see  that  the  wind 
velocities  of  22.292  mph  for  bending  and  44.837  mph  for  torsional  vibrations 
are  most  likely  the  ones  to  be  concerned  about.   The  larger  velocities 
corresponding  to  vibrational  modes  2  and  3,  in  most  areas,  will  not  occur  too 
frequently.  If  22  mph  and  45  mph  winds  are  a  very  common  occurrence  in  the 
location  of  the  structure,  the  designer  may  decide  to  either  modify  the  member 
using  methods  presented  in  this  report  or  choose  a  new  member  for  analysis 
in  order  to  reduce  possible  detrimental  fatigue  damage  due  to  wind  induced 
vibrations. 
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DESIGN  EXAMPLE  NO.  2 

This  example  shall  examine  another  "H  ' shape  section,  this  time  a 
large  built  up  member  with  both  ends  having  fixed  supports.   (o1=a2  ->  <=° 
for  bending  and  torsion) 


Properties  of  member 
And  orientation  in 
the  wind 


V 


E  =  29  x  10-5  ksi 
G  =  1.12  x  104  ksi 


It 

w 


nr 


:f-*i 


Area  =  38.5  in. 

h  =  30.0  in.  (Depth  of  section) 

"  d  =  24.0  in.  (Width  of  flange) 

tf=  0.5  in. 

tw=  0.5  in. 

Ixx=  6238.208  in4 

Iyy=  1152.302  in.4 

:t  I„  =  7390.510  in.4 


E  =  3.229  in. 

Cw=  2.51  x  105  in.6 


The  member's  length  is  90  feet  and  it  carries  a  prevailing  stress  of  30ksi 

The  tensile  load  is  therefore: 

T  =  a   x  A  =  (30ksi)(38.5  in2)  =  1155.0  kips. 

Now  calculate  /P: 

2 
TL  ,_ 

For  bending:   PB=  ~   =  40.31,  so  /P  =  6.35 


yy 


For  torsion:   PT  =  (GC  +  I  T/  )L2  =  41.32,  so  /p  =  6.43 


EC 

EUw 

The  values  of  v^Q  can  now  be  read  from  the  chart  (Fig.  10)  on  page  24. 
Again  using  a  Strouhal  number  of  St=0.12,  (for  "H"  shaped  section),  Z  can 
be  computed  using  Z  =  ^Qt  v/St.  The  results  of  these  efforts  are  sum- 
marized in  the  following  table: 


Mode 

Bending 

Torsion 

n 

H 

ZB 

* 

zT 

1 

5.58 

16.108 

5.60 

16.166 

2 

8.60 

24.826 

8.65 

24.970 

3 

11.65 

33.631 

11.70 

33.775 
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The  velocities  for  the  first  three  modes  of  vibration  can  be  calculated 
using  the  following  relations: 

„   d 


For  bending    V  =  Z 


EI 


n"  Bn  2TrL  J  tf 


¥¥- 


For  torsion    V„=  Z^       ' 


2i\L' 


The  velocities  calculated  are  summarized  in  the  following  table  which  also 

show  the  corresponding  Reynolds  number  calculated  from: 

Re  =  Vnd 
n   n 


v 


Mode 

Bending 

Torsion 

n 

Wind 
Velocity 

V 
MPH(Km/hr) 

Reynolds 

Number 
Re 

Wind 
Velocity 
V 
MPH(Km/hr) 

Reynolds 

Number 

Re 

1 
2 
3 

52.52(84.52) 
124.84(200.91) 
229.102(368.704) 

9.645xl05 
2.292xl06 
4. 208x1 06 

56.39(90.75) 

134.53(216.51) 
246.14(396.13) 

1.04xl06 

2.47xl06 
4.52xl06 

As  previously  explained,  the  designer  at  thin  point  must  decide  on  the 
worthiness  of  the  member  in  light  of  the  above  information  and  act  accord- 
ingly. 
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DESIGN  EXAMPLE  NO.  3 


This  example  will  examine  how  the  theory  is  applied  to  a  steel 
pipe  -  12"  0  standard  weight. 

Properties  of  Member 

O.D.  =  12,750  in. 

I.D.  =  12.000  in. 

Nominal  Diameter  =  12  in. 

Wall  thk.  =  0.375  in. 

Area  =  14.6  in.2 

I  =  279  in.4 

I  =  558  in.4 
C  =  %irr4  =  2035.75  in.4  P 

Cw=  0     ,  Length  =  40  ft. 

E  =29  x  10  ksi. 
G  =  1.12  x  10  ksi. 


Support  Conditions:   Both  ends  simply  supported  (o*  =   a  =  0) 
Loading  Conditions:   Assume  prevailing  tensile  stress  of  35  ksi, 
therefore  the  tensile  load  is: 


T  =  (35  ksi) (14.6  in.  )  =  511  kips. 


Calculations 


1.  Caclulate  P  for  bending  (Note:  P  for  torsion  =  0  because  C   =0) 

2 
P  =  JL 2  =  (511  kips)(40ft)2(144in-/ft2  )=   14.55 

EI      (29  x  10-3ksi)(279  in.4) 
/P  =  3.81 
Using  the  above  value  of  /p~  and  Fig.  1  (p.  15),  the  value  of  /Q  for 
the  first  three  modes  of  vibration  in  bending  can  be  obtained  while 
Table  1  (p.  25  )  provides  the  values  of  -/Q  for  torsion.   Using  the 
values  of  7Q  and  a  Strouhal  number  for  circular  cylinders  of  0.20,  the 
value  of  Z  can  easily  be  calculated  from  the  following  relationship: 

Z  =  $Q  r  ,/sT 

4y— 
The  following  table  lists  the  values  of  /Q  and  Z  for  the  first  three 
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modes  of  vibration  in  this  example, 


MODE 
n 

,    BENDING 

TORSION 

Jq 

Z 

k 

Z 

1 

2 
3 

3.90 
6.75 
9.75 

8.72 
15.09 
21.80 

1.77 
2.51 
3.07 

3.96 
5.61 
6.86 

The  wind  velocities  corresponding  to  the  first  three  modes  of  vibration 
can  be  calculated  using  the  following  relationships: 

For  bending:  ?     d    /  EI 

Vn  =  Zn  2ttL2V  pA 
m 


For  torsion: 


GC 


Vn  =  \   M2J     pi 

11     ^   m  P 


The  Reynolds  numbers  corresponding  to  the  first  three  modes  of 
vibration  can  also  be  computed  once  the  velocities  have  been  determined. 

as  follows:  ,,  , 

V  d 
Ren=  n_ 

v 

-2  '   2 
where  v  =  kinematic  viscosity  of  air  (v  =  2.3  x  10"  in-/    ) 

sec. 

The  following  table  shows  the  velocities  and  Reynolds  numbers  for  this 
example. 


Mode 

Bending 

Torsion 

n 

Wind 
Velocity 
MPH  (km/hr) 

Reynolds 

Number 
Re 

Wind 
Velocity 
MPH  (KM/hr^ 

Reynolds 

Number 
Re 

1 
2 
3 

33.10(53.27) 
99.12(159.52) 
206.88(332.94) 

3.23  x  105 

9.67  x  105 

20.17  x  105 

1.85(2.98) 
3.72(5.99) 
5.56(8.95) 

1.81  x  104 
3.63  x  104 
5.42  x  104 

Note:   For  this  case  torsion  seems  to  be  the  dominating  design  factor 
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In  examining  the  Renolds  number  for  the  first  mode  of  vibration  in 
bending,  it  can  be  seen  that  the  flow  falls  within  the  transition  range 
for  circular  cylinders  as  shown  in  Fig.  17,  P.  4  7.  Therefore,  before 
making  any  final  decisions  concerning  the  member  the  designer  should 
check  for  random  shedding.   This  example  will  be  continued  in  design 
example  No.  5. 
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II.  APPLICATION  OF  HARTLEN-CUKRIE  MODEL  TO 
ELASTIC  BEAM  IN  SUBCRITICAL  FLOW  CONDITION 

From  previously  presented  material,  the  natural  frequency  of 
a  given  beam  subjected  to  wind  forces  can  easily  be  found.   The 
next  important  step  the  designer  must  face  is  the  determination 
of  the  maximum  deflection  the  beam  will  undergo  for  a  certain  mode 
shape  due  to  a  wind  load  and  the  corresponding  stress  produced 
in  the  beam  due  to  this  deflection. 

If  the  deflection  of  the  beam  due  to  the  wind  excitation  is 
negligible,  the  response  can  of  course  be  determined  by  the  conven- 
tional linear  resonance  theory  for  damped  vibrations.   However, 
the  problem  inherent  in  this  approach  is,  if  the  deflection  is 
negligible,  the  associated  flexural  stresses  should  also  be  negli- 
gible.  A  reasonable  approach  must  involve  the  structure-fluid 
interaction.   The  only  well-developed  interaction  model  is  due  to 
Hartlen  and  Currie,  Chapter  II,  Part  B.   Unfortunately,  their 
model  deals  with  a  rigid  beam  while  our  present  interest  is  focused 
on  an  elastic  member  which  is  capable  of  resisting  bending  moment. 
In  order  to  use  their  theory  for  elastic  beams,  a  suitable  hypo- 
thesis must  be  made.   As  shown  in  Fig.  31,  it  is  conceivable  that 
some  segment  of  the  actual  beam  being  analyzed  will  behave  like 
the  rigid  model  with  the  spring  end  conditions  being  provided  by 
the  remaining  portions  of  the  beam. 

Using  the  Hartlen-Currie  Model ,  for  a  given  spring  constant , 
the  deflection  of  the  beam  can  be  found.   Using  the  assumption  that 
some  segment  of  the  beam  under  analysis  responds  as  the  model,  the 


designer  has  merely  to  determine  a  suitable  spring  constant  for  the 
segment.   This  can  be  accomplished  by  the  following  postulation:  the 
natural  frequency  of  the  rigid  beam  is  equal  to  one  of  the  natural 
frequencies  of  the  elastic  beam,  i.e.,  the  natural  frequency  of  the 
elastic  beam  is  equal  to  the  square  root  or  some  constant  (k)  divided 
by  the  mass  density  of  the  beam  ,  .m  , 


*  m 


m 

if  this  natural  frequency  w  is  set  to  equal  a  particular  natural 
frequency  of  the  elastic  beam  under  consideration,  the  constants  k 
can  be  solved  in  terms  of  m  and  w.   It  is  therefore  reasonable  to 
assume  that  this  k  corresponds  to  the  spring  constant  used  in  the 
Hartlen-Currie  model.   Thus  the  designer,  knowing  the  mass  density 
of  the  beam  and  its  natural  frequency/  can  solve  for  the  value  of  k 
and  use  it  in  the  Hartlen-Currie  method  of  analysis  to  determine  the 
deflection. 

It  should  be  noted,  however,  that  the  deflection  determined 
by  the  Hartlen-Currie  model  represents  the  mean  deflection 
of  the  beam.   The  designer  will,  of  course,  be  interested  in  the 
maximum  deflection  for  the  analysis  of  the  beam..   Therefore,  the  de- 
flection resulting  from  the  model  analysis  must  be  modified.  Through 
integration  of  mode  shapes, (see  Table  4),   it  was  found  that  the 
maximum  deflection  of  the  beam  would  be  on  the  order  of  1.5  to  1.7 
times  the  mean  deflection.   Using  this  information  and  a  desire  to 
be  on  the  safe  side,  the  designer  should  consider  the  maximum 
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deflection /W    ,  to  equal  2.0  times  the  mean  deflection/  w,  given 
max 

by  the  Hartlan-Currie  model  analysis  , i.e.  , 


w  _=v  =   2w 

max 


The  determination  of  maximum  deflection  for  known  beam  loading 
conditions  and  wind  forces  is  a  simple  task.   Once  the  deflection 
is  known  the  stress  can  easily  be  found.  When  the  designer  knows 
the  maximum  deflection  of  the  beam  for  a  given  deflection  curve  and 
wind  force,  the  stress  placed  on  the  beam  by  the  deflection  can  be 
determined  using  moment-deflection  relations. 

The  assumptions  made  in  the  above  analytical  method  can  be 
summarized  as  follows :   It  is  first  assumed  that  some  small  segment 
of  the  elastic  beam  being  analyzed  will  act  like  the  rigid  beam  of 
the  Hartlen-Currie  model.   The  remaining  portions  of  the  beam  will 
act  to  provide  the  spring  like  supports  used  for  the  model  analysis. 
Secondly,  it  is  assumed  that  the  damping  and  frequency  of  the  small 
segment  is  similar  to  that  of  the  spring  supported  rigid  beam.   Thus, 
the  natural  frequency  and  mass  density  of  the  elastic  beam  can  be 
used  to  calculate  a  constant ,  k  ,  which  corresponds  to  the  spring 
constants  of  the  Hartlen-Currie  model.   A  third  assumption  made  is 
that,  for  calculating  maximum  deflection, the  designer  use  a  correc- 
tion coefficient  of  2  times  the  mean  deflection  resulting  from  the 
model  analysis.   This  correction  coefficient  may  be  slightly  conser- 
vative, but  it  gives  a  more  accurate  value  of  the  maximum  deflec- 
tion needed  for  design  calculations  than  the  direct  model  analysis 
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results.  The  final  assumption  that  has  to  be  made  is  that  the  wind 
force  acts  uniformly  on  the  beam. 

The  moment  and  the  deflection  in  an  elastically  supported  beam 
can  be  easily  formed  by  the  moment-area  method.  An  elastically  supported 
beam  may  be  considered  as  a  simply  supported  beam  with  an  end  moment 

applied  on  each  end.   It  is  well-known*  that  the  maximum  moment  of 

1   9 
simply-supported  beam  subjected  to  a  uniform  load  is  tt  qL  .   The 

slope  at  the  end  is  3PI  g-  qL2  =  2AEI    •      The  sloPe  due  to  tne  end 

moment  is  M  L__   The  end  moment  MQ  is  yet  of  course  unknown,  but  is 

2EI 
in  turn  related  to  the  slope  at  the  end  to  the  combined  effect 

M0  =  k  dv  =  EK_  dy 
dx    L  dx 

Thus   EI  dy  =  _1_    qL3  -  EI  G_  dy  .   Solving  for  dy  and 
dx   24  2  dx  dx 

substituting  into  the  expression  for  M  ,  we  obtain 

Mo  =      g    °^ 
12(2  +  a) 

The  maximum  deflection  is 

w    =   5  qL4   -  MQ   L2 


max   384EI         8  EI 
The  results  are  tabulated  in  Table  3,  and  plotted  in  Figure  14 
In  the  application,  the  following  data  are  usually  given: 

a  ,   d,   Lj   in,   and  y  . 
It  is  desired  to  compute  the  flexural  stress  f,  . 


*For  simplicity  we  restricted  our  attention  to  a  common  case  where  the 
beam  is  symmetric  and  prismatic  and  is  uniformly  loaded. 
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Table  3.  Deflection  Moment  Relationship 


End  Fixity 
Coefficient  ,<? 

Maximum  Moment* 
Coefficient^ 

Maximum  Deflection 
Coefficient ,b 

a 
b 

0 

0.125 

0.01302 

9.600 

0.1 

0.121 

0.01253 

9.657 

0.5 

0.1083 

0.01094 

9.899 

1 

0.0972 

0.00955 

10.178 

5 

0.0655 

0.00557 

11.759 

6 

0.0625 

0.00521 

11.996 

10 

0.0695 

0.00432 

16.088 

20 

0.0758 

0.00352 

21.173 

50 

0.0801 

0.00302 

26.523 

100 

0.0817 

0.00282 

28.972 

200 

0.0825 

0.00272 

30.330 

1000 

0.0832 

0.00262 

31.756 

10000 

0.0833 

0.00261 

31.916 

1000000 

0.0833 

0.00260 

32.039 

00 

0.0833 

0.00260 

32.039 

Use  value  of 


12  (2  +o   ) 


or   -  - 


12(2  +a  ) 


whichever  is  larger. 
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The  procedure  is  as  follows :   For  given  Y  ,   one  can  compute 

YAPm  _ 

the  nondimensional  damping —  ,from  which  one  can  determine  w/d, 

pad2  ' 

according  to  Hartlen-Currie  model.   From  elementary  beam  theory, 

bqL4  2 

w    = — £- —  ,  M  =  aql/ 

max    EI 

Eliminating  q  between  them, 


M  =  f  E  w   — 
b   max^2 

By  flexural  formula, 

W  jo 

„  M  c    a     max    d  2 

f  = =   -   E  (  -  ) 

b   I      b     2d      L 

Noting  that  6/2=5  ,    we  finally  obtain      w   n  =  w 

max' 

a    w    d  2 
fh=  g  E  -   (  -  ) 
b   b    d    L 

where  a/b  corresponding  to  a  given  a  can  be  obtained  from  Fig.  14 


DESIGN  EXAMPLE  NO. 4:  Using  the  cross  section  of  example  1  (p. 116) ,  the  mean 
deflection  and  flexural  stress  for  a  member  in  the  subcritical  regime  can  be 
calculated. 

Given  that  the  length  (L)  =  100  ft,  end  fixity  a-,  =02=10,  structural  damping  co- 

-7      2    4 
efficient  (y)  =  0.0056,  and  mass  density  of  air  (p  )  =  1.123x10   lb-sec  /in 

a.  Mean  deflection  can  be  found  from  Fig.  11,  p.  26 

2 
yApm/pad  =  5„33  ;  therefore,  w  =0.09  which  means  w  =  0.0074  inch. 

~d~~ 

b.  Flexural  stress  can  be  found  as  follows:  From  Fig.  14,  p.3o.  a/b  =  16 

f ,  =  a  E  w   (d  )2  =  16  x  29  x  106  x  0.09  x  ,12.217   ,2  =  4328  psi 
&   d    L  100  x  12 
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rms  w 


III.    APPLICATION  OF  RANDffi    VIBRATION  THEORY 
TO   SUPERCRITICAL    FLOW  CONDITION 

Equation    (7)    in  p.    109   can  be  wriiten  as   follows: 
1  p.d 


L 
x=- 

2 


le/rr3 


CT    -a-   V2    (L.    +   L_    +   L      + ) 

L   n    A  1  2  3 


0.075    lb/cf 


487   lb/cf 


m 


(for   steel) 


L      rL 


L.    = 

l 


wi(o)  b  widx   /sJTstryst. 


f2   'L  w2  dX 


i  o 


6j/2tt 


where  Yt  ^s  substituted  by  5./2tt. 

For  design  purposes,  we  shall  assume  5-  =  6  =  constant. 

Now,  for  symmetric  beams,  L„  =  0 ,    and  the  product  of 


w.(|)  'L 

-^-   and  -2- 


w±  dX 


(where  K 


K. 


K. 


i  =  x/i  Wi 


dX  ) 


can  be  found  in  Table  4.   It  is  observed  that  each  is  virtually 
independent  of  a  and  P  .   Thus,  the  products  are  slightly  less  than 
1.4  and  0.5  for  i  =  1  and  i  =  3,  respectively. 
We  observe  that, 


Sf(St1)St    1.4 

in  which,  the  values  of  S<?,  can  be  obtained  from  Figs.  12  and  13.   It 
is  noted  that  S£  (St3)/Sf  (St-, )  and  St  /St,  are  approximately  0.2  and  3, 
respectively.   Thus  for  f^/f  <  0.2,  —  <  1%. 
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Table  4  Mean  and  RMS  Deflection 
First  Mode  (n=l) 


°1  =°2 

Jv 

4?])/ki 

Vi>/  Kl 

o.o 

.90 

1.4 

0 

7.7746 

.90 

1.4 

10.9545 

.90 

1.4 

0.0 

.88 

1.5 

5 

7.7746 

.89 

1.5 

10.9545 

.89 

1.4 

0-0 

.87 

1.5 

10 

7.7746 

.88 

1.5 

10.9545 

.88 

1.5 

0.0 

.86 

1.5 

20 

7.7746 

.87 

1.5 

10.9545 

.87 

1.5 

0.0 

.84 

1.6 

100 

7.7746 

.85 

1.5 

10.9545 

.86 

1.5 
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Table  4  (Con't)  Mean  and  RMS  Deflection 


Second 

Mode  (n=2) 

ol  =  o2- 

J* 

3^2VK2 

w2(|)/K2 

0.0 

0 

0 

0 

7.7746 

0 

0 

10.954S 

0 

0 

0.0 

0 

0 

5 

7.7746 

0 

0 

10.9545 

0 

0 

0.0 

0 

0 

10 

7.7746 

0 

0 

10.9545 

0 

0 

0.0 

0 

0 

20 

7.7746 

0 

0 

10.9545 

0 

0 

0.0 

0 

0 

100 

7.7746 

0 

0 

10.9545 

0 

0 
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Table  4  (Con't)  Mean  and  RMS  Deflection 
Third  Mode  (n=3) 


Oi=a2 

Vp 

^K3 

w3(^)/K3 

0.0 

.30 

-  1.4 

0 

7.7746 

.30 

-  1.4 

10.9545 

.30 

-  1.4 

0.0 

.33 

-  1.4 

5 

7.7746 

.32 

-  1.4 

10.9545 

.31 

-  1,4 

0.0 

.35 

-  1.4 

10 

7.7746 

.33 

-1.4 

10.9545 

.32 

-  1.4 

0.0 

.36 

-  1.4 

20 

7.7746 

.34 

-  1.4 

10.9545 

.33 

-  1.4 

0.0 

.37 

-  1.4 

100 

7.7746 

.34 

-  1.4 

10.9545 

.33 

-  1.4 

L 
Where       ^  \  =  /  w  dX 


n'    o   n 


L  2  h 
K   ■={   f     w  dX) 

o   n 


n  =  1,  2,  3. 
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Then  L3  is  less  than  1%  of  Lj  and  L^  should  represent  at  least  95%  of 

00 

£  L^ .      Under  this  circumstance, 

i  =   1,3 


rms  w     =   2.445  x  10  6     V2 


w:    (2)    I(wi)  /  SfCSt^St! 


Aff       L      Kl  h 


cSj/tt 


where  I(w-j_)   denotes       •>      w-,dx, 


=   3,1  x  10        V2 

Af1 


6i/tt 
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DESIGN  EXAMPLE  NO. 5 

This  example  is  a  continuation  of  example  No.  3.  Here  the  member 
shall  be  checked  for  random  vibration.  Example  No.  3  showed  that  the 
Reynolds  number  for  the  first  mode  of  vibration  in  bending  fell  within 
the  transitional  range  of  1.0  x  105  to  4.0  x  10  as  indicated  in  available 
literature  (see  Fig.  17  p.  47  and  Fig  19  p.  50).  Both  references  deal 
with  circular  cross  sections*,  and  so  an  analysis  shall  be  made  on  the 
member  of  example  No.  3. 

Beam  properties  are  listed  on  p.  123.  12" $  standard  weight  pipe  . 
Both  ends  are  simply  supported  (o^  =a2  =  0) 

Length  =  40  Feet  Tensile  Stress  =  35ksi 

Independent  Lift  Coefficient  Cy  =  1.2 
Logarithmic  Decrement  5 1=<$2  =  0.02 

From  example  No.  3  the  velocity  corresponding  to  the  first  mode  of 
vibration  in  bending  is  33.10  mph  which  equals  582.57  in/sec.  The  actual 
frequencies  of  vibration  for  the  first  three  modes  of  vibration  can  easily 
be  calculated  knowing  the  values  of  «/Q  from  the  following  relationship: 

Q4  EI  "* 
m  A  L4 
Where     P  =  Frequency  (Radians/Sec)     n  =  1 .  2 ,  3 

4 

n  =  v^q 
Pro  =  Mass  Density 


*Note:  The  Reynolds  numbers  for  HP  sections  of  examples  1  §  2  also  fall 
above  the  1.0  x  10^  zone  limit  however,  no  literature  is  available  to 
indicate  that  cross  sections  other  than  circular  ones  reach  supercritical 
vibrations.   If  it  is  found  that  structural  shapes  do  infact  reach  a 
supercritical  range  then  the  method  outlined  in  this  example  should  be 
used  to  check  flexural  stress  due  to  random  vibration. 
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Thus  for  this  example 

P  =  57.42  Rad/Sec 

P2  =  172.02  Rad/Sec 

P3  =  358.90  Rad/Sec 

These  frequencies  can  now  be  used  to  find  the  normalized  power 
spectrum  S£(Stn)  which  is  a  function  of  the  resonant  Strouhal  number* 

Stn  =  


P 


2ttV1 

Therefore  St..  $  St,  can  be  used  to  find  S_  (Sti)  and  S  (St,)  on  Fig.  12 
27  i     o  f    x       f   ° 

Stx  =  0.200  SfCSt1)  =  1.3 

St3  =  1.250  Sf(St3)  =  0.1 


From  Table  4 


w,(L/2) 
^  =  -  1.4 


L 


w 

/   -1  dx  =  0.90 

0   K1 


L 
0 


w3 
/   -  dx  =  0.30 

K 


* 


Note:  Since  the  Reynolds  number  is  within  the  transitional  range  the 
vibration  could  become  harmonic,  that  is,  the  wind  velocity  for  the 
1st  mode  excites  vibration  in  the  third  mode.   Therefore  for  random 
vibrations  Strouhal  numbers  must  be  calculated  for  the  1st  and  3rd  modes 
using  the  velocity  of  the  1st  mode. 
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With  the  above  information  the  root  mean  square  deflections  of  the 

member  at  irddspan  can  be  calculated  using  the  formula  presented  on  p.  1 35 

as  follows: 

For  n  =  1 : 

(582.57  in/sec)2 
rms  w(        =  (2.445  x  10~6)(1.2) 


x  =  L/2  /57.42  rad/sec\2 

(14.6  in2)} 


/57.42  rad/secV 
(  2,  ) 


(12.75   in)(1.4)(0.9) 


=  0.084   in 


For  n  =   3 


rms  wi  =    (2.445  x   10"°)  (1.2) 


(1.3) (0.20)   ' 
(0.02/  tO 

/ 
(582.57  in/sec)2 


x  =   L/2  ,358.90  rad/sec%z 

(14.6  in2)  ( -1 

^      2  TT 

|(0.1)(1.25)  s 
(12.75  in)(-1.4)(0.30) 


(0.02J/7T 


=  -  4.96  x  10~4  in 


So  from  the  above  it  can  be  seen  that  the  deflection  of  the  member 
due  to  resonant  vibration  action  tends  to  reduce  the  principal  deflection 
due  to  the  natural  vibration  at  the  first  mode  by  0.6%,  a  negligible 
amount  for  this  example.  The  actual  midspan  deflection  is  then: 
w  =  (8.4  x  10-2  in)  -  (4.96  x  10"4  in)  =  8.35  x  10"2  in 
The  bending  stress  (f^)  for  this  deflection  can  easily  be  calculated  as 


f ,  =  r   E  -  (;h2       where:  f  can  be  read  from  Table  3  p.  134  or 

found  in  Figure  14  p.  31 


So:  fh   =  (9.6) (29  x  106  lb/in"2) 


(For  Oi=  a2=  0  a/b  =  9.6) 
(8.35  x  10~2  in)  /     (12.75  in)   N? 
(12.75  in)   \(40  ft)  (12  in/ ft)/ 


i 


=  1.286  x  103  lb/in2  =  1.286  ksi 
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It  is  interesting  to  compute  f,  by  using  the  relation 


2 

9  w 
Bending  Moment  (M)  =  -EI  — 

3  x: 


Since  w  =  Sin  x  it  follows  that 


L 

X=  2 


2 

9  w       T 
9  x2  I  x  =  - 

2 


w 


and 


=  EI 


=  (29  x  106  lb/in2) (279  in4) 


=  2.894  x  10  in  lbs 


(40  ft) (12  in/ ft) 


(8.35  x  10"2  in) 


Mc 


and  f,  = 
b 


(2.894  x  KT  in  lbs)  (12.75  in) 
(279  in4) (2) 


=  6.613  x  10  psi 


?  3 

At  peak  frequency  f,  =  (2)  (6.613  x  10  psi)  =  1.3  x  10  psi 

which  checks  with  the  previous  calculation  of  f,  . 

With  this  knowledge  the  designer  now  has  a  complete  picture  of 

how  the  chosen  section  will  react  under  existing  winds.  A  sound 

decision  can  therefore  be  made  as  to  whether  to  use  the  given  member 

as  it  is,  modify  it  to  alter  detrimental  vortex  shedding  patterns,  or 

choose  a  new  section  for  analysis. 
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PART  D 
BIBLIOGRAPHY 


The  main  approach  to  search  for  relevant  literature  on 
flow-induced  vibrations  is  as  follows: 

1.  Personal  Collections 

Dr.  Nicholas  Basdekas,  is  a  specialist  in  this  area,  and 
has  been  interested  in  this  subject  for  two  decades.   In  his 
professional  activities,  he  has  maintained  contact  with  a 
large  number  of  researchers  in  this  field.   He  has  kindly 
loaned  us  his  personal  collection  of  thirty  articles  in  this 
area. 

Dr.  Chris  Yiu  of  Pavlo  Engineering  Co.  of  New  York,  is 
currently  the  Chairman  of  Subcommittee  on  Cable  Stayed 
Bridges  in  the  ASCE  Structural  Division.   His  Committee  is 
the  clearinghouse  of  information  on  the  cable  stayed  bridges 
and  has  recently  produced  a  draft  copy  of  bibliography  on 
Cable  Stayed  Steel  Bridges.   In  this  collection,  it  contains 
13  papers  on  aerodynamic  behavior  of  cable  stayed  bridges  and 
12  papers  on  behavior  of  cables. 

2 .  Library  Search 

Chi  Associates  has  access  to  the  facility  of  computerized 
literature  search.   The  following  sources  were  scanned  for 
the  period  covered  in  parathesis: 

*  National  Technical  Information  Service  (1964-1975) 

*  Engineering  Index  (1970-1975) 

*  Defense  Document  Center  -  Unclassified  Reports  (1966- 
1975) 

*  INSPEC-MECH  ENGR  (197  3-1975) 
The  following  key  words  were  used: 

*  Strouhal 

*  Strumming 

*  Vortex  Shedding 

*  Vortex  Street 

*  von  Karman  Street 


*  Flow  Induced  Vibrations 
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The  number  of  hits  scored  in  each  key  word  for  each  source 
was  recorded, and  all  hits  were  identified  in  terms  of  the 
titles  of  the  articles  or  report  and  the  abstracts,  if  avail- 
able.  Through  the  review  of  the  titles  and  abstracts,  a 
decision  was  made  which  of  these  articles  or  reports  should 
be  acquired  for  closer  examination.   All  articles  processed 
this  way  are  listed  in  alphabetical  order  of  the  first  au- 
thor's names. 
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Novo!  fteaeej-efo  leibstrateiry 

Washington,  D.C.  20390 

13  August  1975 


Tel.  202  767  2220 

8404: 162 :RHV:bgs 


Dr.  Michael  Chi 

2721  24th  Street  North 

Arlington,  VA  22207 

;ar  Dr.  Chi: 

This  is  reply  to  your  request  for  information  about  the 
Strouhal  number  for  structural  shapes.   A  search  of  tb 
Shock  and  Vibration  Information  Center's  data  base 
performed  using  the  following  key  words  or  the  combination 
of  key  words : 

Vortex  Shedding 

Vortex  Noise 

Vortex  Induced  Excitation 

Structural  Forms 

Structural  Members 

No  references  to  reports  that  would  contain  values  of  the  Strouhal 
number  for  structural  shapes  were  found  in  this  search .   The  Shock 
and  Vibration  Information  Center's  data  base  contains  references 
to  many  technical  papers  and  unclassified  unlimited  technical 
reports  that  have  been  published  since  1969. 

We  did  find  a  reference  to  the  range  of  Strouhal  numbers  in  a 
paper  by  Scanlan  and  Wardlaw  in  the  following  reference: 

Scanian,  R.H„  and  Wardlaw,  R4L.,  Reduction  of  31   -Induced 
Structural  Vibrations,  Proceedings  of  Colloquium  on 
Isolation  of  Mechanical  Vibration,  Impact,  and  Noise, 
(AMD-1),  presented  at  the  Design  Engineering  Technical 
Conference,  Cincinnati,  Oh,  September  1973,  p  41. 

Very  truly  yours. 
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Principal  Review  Articles : 

Morris,  A.  W. ,  "A  Review  on  Vortex  Streets,  Periodic  Wakes 
and  Induced  Vibration  Phenomena,"   J.  of  Basic  Engineering, 
Trans.  ASME,  (June,  1964) ,  pp.  185-196. 

Berger,  E.,  and  Wille,  R. ,  "Periodic  Flow  Phenomena," 
Annual  Review  of  Fluid  Mechanics,  Vol.  4.,  1972. 

"Wind  Forces  on  Structures,"   Final  Report  of  Task  Committee 
on  Wind  Forces,  Trans.  ASCE,  Vol.  126,  Part  II,  1961,  pp.  1124- 
1197. 

Diggs,  J.  S.,  "A  Survey  of  Vortex  Shedding  from  Circular 
Cylinders  with  Applications  Toward  Towed  Arrays,"   Report  by 
MAR,  Inc.  to  NSRDC,  Cardarock,  Maryland,  1974. 

Wille,  R. ,  "Karman  Vortex  Streets,"   Advances  in  Applied 
Mechanics  VI ,  Academic  Press,  New  York,  1960. 

Wille,  R. ,  "On  Steady  Flows  and  Transient  Motions,"   Progress 
in  Aeronautical  Sciences,  Vol.  7,  Pergamon  Press,  New  York, 
1966. 

Morkovin,  M.  V.,  "Flow  Around  Fully  Separated  Flows,"  ASME 
Symposium  on  Fully  Separated  Flows,  1964. 

Annotated  Bibliography: 

Mulcahy,  T.  M. ,  and  Chen,  S.  S.,  Annotated  Bibliography  on 
Flow  Induced  Vibrations,   Argonne  National  Lab.  Tech.  Memo. 
ANL-CT-74-05. 

Symposia  and  Books : 

The  Modern  Design  of  Wind-Sensitive  Structures,  CIRIA  Publication, 
Proc.  of  the  Seminar  held  on  June  18,  1970,  at  ICE  London. 

Proceedings  of  International  Symposium  on  Vibrational  Problems 
in  Industry,  1971. 

Sachs,  P.,  Wind  Forces  in  Engineering,  Pergamon  Press,  New 
York,  1972. 

Flow-Induced  Structural  Vibrations,  IUTAM-IAHR  Symposium, 
Karlsruhe,  E.  Naudascher  (Ed),  Springer-Verlag,  New  York,  1972. 

Proc.  of  the  Third  International  Conference  on  Wind  Effects 
on  Buildings  and  Structures,  Tokyo,  Saidon  Co.,  Tokyo,  1971. 
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ASME  Symposium  on  Fully  Separated  Flows,  1964. 

Mair,  W.  A.,  and  Maull,  D.  J.,  "Bluff  Bodies  and  Vortex 
Shedding  -  A  Report  on  Euromech  77,  J.  Fluid  Mech. ,  Vol.  45, 
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PART  E  (APPENDIX  B)    COMPLETE  SET  OF  GRAPHS 
FOR  NATURAL  FREQUENCY  COEFFICIENTS  AT  DIFFERENT  LEVELS 

OF  AXIAL  TENSION 
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